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SUMMARY 

NASCAP is a computer program that comprehensively 
analyzes problems of spacecraft charging. Using a fully 
three-dimensional approach, it can accurately predict 
spacecraft potentials under a variety of conditions. Under 
Contract NAS3-2176:2, several changes were made to improve 
NASCAP, and a new code, NASCAP /LEO , was developed. In addi- 
tion, detailed studies of several spacecraft-environmental 
interactions and of the SCATHA spacecraft were performed 
(see companion reports. References 1 and 2) . 

The changes made to NASCAP this year make it run 
faster - execution time has been cut by more than half. It 
is now more accurate - photosheath approximations and many 
other features have been added. It is easier to use — a 
new user unfamiliar with the host computer can learn NASCAP 
in a matter of days. Auxiliary codes have also been de- 
veloped to improve NASCAP flexibility and reduce problem 
turnaround time. 

A NASCAP/LEO program has been written. This handles 
situations of relatively short Debye length encountered by 
large space structures or by any satellite in low earth 
orbit (LEO) . This preliminary model has already been applied 
to two important high voltage interaction problems. 
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1. INTRODUCTION 


This is the Pinal Report on Tasks 5 and 6 of Con- 
tract NAS3-21762, "Additional Application of the NASCAP Com- 
puter Code", The work was performed by Systems « Science and 
Software between 26 March 1979 and 9 OctoJ^er 1980. 

This report includes much of the information from 
monthly reports published throughout the contract, as well 
as some new items. There are two other volumes of the final 
report. One contains environmental interaction studies of 
three space systems. Tasks 1, 2, and 3; the other describes 
work done on Task 4 — SCATHA spacecraft study. 

NASCAP ( NASA Charging Analyzer Program) is a powerful 
computer program for studying spacecraft charging. It is 
intended as a tool for spacecraft designers to help them 
avoid the problems of electrical charging in space. It has 
also been used by scientists associated with the SCATHA satel- 
lite to interpret the results of on-board experiments. 

Since there is now a group of researchers using NASCAP 
on a day-to-day basis, there is now a need for code extensions, 
improvements, and miscellaneous modifications. The changes we 
have made to NASCAP this year fall into five categories: phy- 

sical model extensions, usability, code review, auxiliary 
codes, and NASCAP/LEO. 

One way to improve the accuracy of a model code is to 
expand the range of physical situations it can handle. Chap- 
ter 2 reports on the changes to the physics of NASCAP. This 
includes new treatments of the photosheath, and improvements 
in detector and emitter modeling. 

Chapter 3 reports on changes that make NASCAP easier 
to use. These changes, along with the new auxiliary codes 
of Chapter 5, make it much easier for a new user to learn 
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NASCAP. More than that, they free the experienced user to 
think more about the actual problem and less about the de- 
tails of code mechanics. The changes involve simplified 

I 

input, more informative output, and a few new features, | 

The most visible of the changes allows all NASCAP input to | 

be in free format - input parameters are no longer card- | 

column sensitive. | 

I 

chapter 4 gives the details of a comprehensive code j 

review. We have made NASCAP faster and cleaner. A parti- I 

I 

cular problem that used to take eight hours of computer time 
now takes two hours. We also put in comments to make the 
source code clearer and made minor improvements. NASCAP 
now can handle up to 1250 surface cells and can handle 
certain shadowing cases that formerly were prohibited. 

Chapter 5 is titled "Auxiliary Codes". These are 
codes to be used in conjunction with NASCAP. In general, 
each auxiliary code takes a single isolated function of 
NASCAP, and performs it easier, faster, and cheaper. 

MATCHG, used before NASCAP is run, allows the user to try 
various combinations of material properties and environments, 
interactively. TERMTALK is an interactive post-processor 
that culls surface cell potential histories from NASCAP 
files. Three more codes — OBJCHECK, P0TPI*T, and FILES - 
were designed for new users in a workshop situation. These 
are not interactive, but provide turnaround an order of 
magnitude faster than otherwise possible. 

The most important change in the NASCAP model is the 
development of an alternate code NASCAP/LEO. This is 
described in Chapter NASCAP analyzes standard sized 
objects in deep space or in geosynchronous orbit. For very 
large objects, or those in low earth orbit (LEO), the NASCAP 
physical models are inadequate. NASCAP/LEO is designed to 
handle these situations. It is a brand new code, and 
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presently has none of the sophistications or convenient user 
features that have been built into NASCAP over the years. 
NASCAP/LEO calculates an equilibrium solution to high voltage- 
dense plasma interaction problems. It handles objects that 
are built from cubes. The code has already been used to study 
ion focusing^ and to calculate the effects of solar voltage 
array patterns on plasma power losses. 

NASCAP/LEO is not ready for general distribution, 

NASCAP is currently available through COSMIC. 


Z. PHYSICAL MODEL EXTENSIONS 


This chapter describes several modifications or exten- 
sions to the physical models employed in NASCAP to calculate 
the charging response of an object in space. The changes 
described below represent fine tuning of existing physical 
models rather than fundamental revisions. 

Enhancements in the treatment of space charge around 
the NASCAP objects due to photoelectrons or the ambient 
plasma are described in Sections 2.1 and 2.3. Updated de- 
fault material parameters are given in Section 2.2. Mis- 
cellaneous modifications to the ROTATE^ DETECTOR, and EMITTER 
routines are described in Sections 2.4 through 2.6. Sec- 
tion 2.7 documents the finite element interpolants used in 
cells containing thin booms. 

2.1 PHOTOSHEATH 

NASCAP has incorporated two new options for calculating 
effects of the photosheath. They are: effective photosheath 

conductivity and self-consistent photosheath. 

The effective photosheath conductivity causes a change 
in surface conductivity. Xt depends on rate of photoemission, 
external electric field, and surface cell voltage. 

The self-consistent photosheath is a feature to be 
used only in very specialized applications. It requires 
certain preconditions and is expensive to use. 

2.1.1 Effective Photosheath Conductivity 

A formulation for effective surface conductivity by 
low energy electrons has been incorporated into NASCAP. The 
option may be invoked by specifying the option 'EFFCON ON' 
and turned off by specifying 'EFFCON OFF'. (The default 
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value IS 'EFFCON OFF’.) The photosheath conductivity is added 
to the intrinsic (material property) surface conductivity at 
each cycle. It is not operative for those cells for which 
surface conductivity is to be ignored. 


The value of the effective photosheath conductivity is 

2 

dependent on J (current of low energy emitted electrons, A/m ) , 
E (external electric field, volts/meter) , V (cell potential), 

M M .3 V ^ W .. A.3 a «... M W ^nW A ...... J J mm J ^ mm mm t .. 1... .. 1 \ m mm 


and X, the code mesh size. The conductivity, a (ohms ) is 
given by 


For E < 0 


- I J x2 


For E > 0 and V > 0 


o = (J/E^) 


( 4 -( 4 - 


2V) 


For E > 0 and V < 0 


a - (J/E^) ^4-(4-2EX) | 


Additionally, a is forced into the range 


lO"^ > a > lO"^® 


2.1.2 Self-consistent Photosheath 


One aspect of spacecraft charging that has not been 
thoroughly investigated using NASCAP has been the structure 
of the space charge barrier formed by photoelectrons. This 
barrier, while small in voltage compared to substorm induced 
charging, is of interest both theoretically and experimentally. 
Experiments , such as SC-10 and SC-2 can make inferences about 





the sheath structure. The extent and magnitude of the sheath 
can be Important in unfolding ambient electric field data. 

Previous NASCAP calculations did not take the space 
charge of the photosheath into account because the voltage 
perturbations of the sheath are extremely small compared to 
the kilovolt/meter fields set up by differential or space- 
craft charging during substorms. However, for the case of 
the unfolding of SC-10 during quiescent conditions, the 
space charge sheath is the dominant source of field differen- 
tials between the two halves of the dipole. This experiment 
has provided electric field strength as a function of dis- 
tance in front of SCATHA. 

In order to predict the space charge barrier, a slight 
modification of the NASCAP explicit photosheath treatment was 
necessary. The 'SHEATH' option in NASCAP predicts photocharge 
densities, but does not use them for potential calculations. 
The space charge calculations run for SCATHA were made using 
the calculated space charge, and the potentials and charge 
densities were iterated on until self-consistency was obtained 
between the potentials and the charge densities obtained from 
particle tracking. This procedure converged quite rapidly. 
Another improvement in the 'SHEATH' routine, that is, emission 
of photoelectrons at several angles rather than just normal to 
the surface, was made to calculate more accurate charge den- 
sities. 

These modified SHEATH routines are not designed for 
general use with NASCAP for several reasons. First, the 
self-consistent routines are presently implemented only for 
fixed spacecraft ground potentials, since the net photo- 
sheath currents are not available to the LONGTIMESTEP features 
of the code. Secondly, to track particles and iterate on the 
potentials would make this procedure prohibitive to use on a 
production basis because of computer time requirements. 
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However, for those interested in scientific investigations of 
the low energy sheath structure the keyword in the RDOPT file 
is 'SHEATH SELF CONSISTENT'. A number can be included on the 
card after these keywords to specify the charge density re- 
laxation parameter, a; 

f ' “"new + 

A value of unity for a would correspond to explicit itera- 
tion. The default value is 

“default " 0*5 

With the default value of the relaxation parameter a typical 
run converged to within ±5 percent in the charge density 
after five iterations. 

For an example of the use of this feature, see 
Schnuelle, et al. 

2.2 MATERIAL PROPERTIES 

In NASCAP, surface material properties are in a state 
of frequent revision. Earlier this year, we conducted a com- 
prehensive reworking of material property parameters. These 
were incorporated into the default values. We are again re- 
structuring our material parameters. 

One of the major changes during the contract year was 
a revision of subroutine PROSEC to give reduced proton second- 
ary yields for protons incident at 1 keV and below. 

Another change was to incorporate emission coefficients 
for protons on metals from Reference 3. Table 2.1 shows the 
new values of NASCAP properties li and 12 for best fit to data 
between 10 and 50 keV. 
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TABLE 2.1. SUGGESTED NASCAP MATERIAL PARAMETERS 
FOR SIX METALS 


Metal 

Property 11 

Property , 

Li 

.335 

91. 

A1 

.244 

230. 

Cr 

.247 

250. 

CU 

.257 

320. 

Ag 

.49 

123. 

Au 

.413 

135. 

(average) 

.33 

191. 


2.3 SCREENING LENGTH 

NASCAP now correctly takes into account the interaction 
between screening length and boundary conditions. It also 
performs a capacitance rescaling when a sudden change of 
screening length takes place. (This procedure is less exact 
than a fresh call to CAPACI, but is far faster and more flex- 
ible.) An interactive program is available to insert neces- 
sary information into files created prior to September 1979. 

2.4 ROTATING MAGNETIC FIELD 

ROTATE now causes any constant magnetic field to be 
rotated along with the sun direction. Any magnetic dipoles 
are not rotated, since they are assumed to be generated by 
the object, and hence would be locked to the object's refer- 
ence frame. ROTATE assumes that the initial, unrotated mag- 
netic field has been specified using the BFIELD keyword in 
a preceding RDOPT. 
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2.5 DETECTOR UPDATE 

The NASCAP detector routines have been modified to 

include a number of new features. This section is intended 

to augment the detector routine documentation contained in 

Katz, I. , et 

Summary of New Features ; 

1. Particle tracking is now performed in all NG computa- 
tional grids of NASCAP. 

2. Particle tracking beyond the highest computational grid 
may. be performed using a 1/R monopole potential. 

3. New I/O buffering scheme permits virtually unlimited 
numbers of particle trajectory line segments to be 
plotted . 

4. The correct potential interpolation function is now em- 
ployed to obtain the electric field for particles which 
pass into boom cells. 

5. Routines have been implemented which correctly determine 
when a trajectory terminates on a boom. (Prior to this 
booms were transparent to NASCAP particle detectors.) 

6. Particle detectors may now be placed on booms. 

7. The energy aperture of NASCAP particle detectors has 
been modified (by popular request) to provide constant 
% energy AE/E width in preference to the old fixed AE 
for all energies. 

8. In addition to the standard energy flux plot produced 
by NASCAP detectors a point-by-point breakdown of the 
flux may (optionally) be printed. 


JE '> 
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9. Optional detailed printout of each discrete particle 
trajectory calculation is available. 

10. Optionally omit flux and trajectory plot headings for 
quick on-line plotting display terminals. 

11. Optional user modification of particle charge and mass. 
Permits tracking heavy ions, for example. 

12. Particle trajectory 2-D projection views may optionally 
be plotted more than once, each time using a different 
maximum grid boundary. This permits the fine details 
of particles returning to the satellite to be viewed 
while also displaying the details of particles circling 
in large radii due to the presence of a magnetic field, 
for example. 

13. All detector keyword input may now be made using free 
format. (The user should still pay attention to key- 
word variable types since a decimal point or exponent 
specified for an integer variable will result in 
errors.) 

Placement of Particle Detectors on Boom Cells ; 

In addition to satellite object surface cells NASCAP 
now permits the user to place a particle detector on a boom 
cell. Placement of a detector on a boom cell is done in a 
way very similar to that used for surface cells. The NASCAP 
user must first have performed an OBJDEP run on the boom- 
containing object which he intends to use. Then referring to 
the boom surface cell list printed by the OBJDEP run he can 
choose the cell number of the boom cell on which he desires 
to place a detector. This cell number is then placed on the 
ICELL keyword card in the detector keyword file. The detector 
will, be located at the axial center of the boom surface cell 
specified by ICELL. In addition, the user must specify the 
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direction of the vector normal to the surface of the boom 
cell. This vector is restricted to lie in that plane which 

a 

is perpendicular to the boom axis and passes through the boom 
surface cell axial center. The vector is specified in terms 
of the angle ? (ZETA) as illustrated in the following diagram 
for a boom along the z-axis. The definition of C for booms 
lying along either of the other two coordinate axis may be 
obtained by cyclic permutation of axis labels in the figure. 
(Note that the coordinate axis are those of the Satellite 
Coordinate System. ) 


Boom Cell 
Section 


Detector on surface 
of boom. Z -coordinate 
of detector location 
is axial center of 
boom cell. 


Normal Vector ' 
to Cell Surface 


Figure 2.1. Particle detector on boom cell. 




The boom cell normal vector defined by the selection 
of ICELL and ZETA plays the same role as the surface cell 
normal vector which defines the Surface Cell Coordinate System 
as described on p. 145-146 of Reference 4, the only dif- 
ference being that the cell normal vector may point in any 
direction within one of the three coordinate planes rather 
than being restricted to one of the 26 mesh symmetry direc- 
tions listed in Table 8.1, p. 147. 

Description of New Detector Parameter Options by Keyword ; 
Detector Independent Parameters 

LABABV Inclusion of this card in a detector keyword defini- 
tion sequence will cause the detector energy flux 
and particle trajectory plot heading labels to be 
suppressed. This feature is useful to the user 
who has access to a Tektronix graphics display 
terminal and desires to speed up his on-line plot- 
ting process. If multiple detectors have been de- 
fined, then this option is not recommended since 
label omission in such cases is apt to be confusing. 

PRBMAT Inclusion of this card causes the inverse matrices 
relating grid potentials to interpolating function 
coefficients in boom cells to be printed. (Nc . 
intended for user applications.) 

PRTRAJ Modulo number for printout of particle trajectory 
steps. (This variable must be an INTEGER even 
though the variable type indicates a real para- 
meter.) If this card is included then for every 
th 

PRTRAJ point of each particle trajectory a line 
will be printed containing the trajectory step 
number (ISTP) , index of NASCAP computational grid 
(IR) , index of boom whose cells the particle is 
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currently in (if it is in a boom cell) (IB) , co- 
ordinates of particle (X) in inner mesh units, 
particle velocity (V) in inner mesh units/timestep, 
electric field (E) in volts/inner mesh unit, and 
magnetic field (B) in gauss. Also printed is the 
electrostatic potential (VBOOMC) which is only 
valid if the particle is within a boom cell 
(IB > 0) . If IB * 0 then the value given for 
VBOOMC is not relevant. 

! The user is cautioned against careless use of \ 
the PRTRAJ option since there is a potential fori 
a tremendous amount of printed output - namely a\ 
maximum of about N*NSTP lines where N is the | 
number of independent variable points! I 

PRPLUX Inclusion of this card will result in a tabulated 
listing of the detector energy flux values in 
addition to the standard plot. (If the PRTRAJ 
option is concurrently active then the listing 
will be interspersed with particle pushing trajec- 
tory information - a possibly undesirable condition.) 

Detector Dependent Parameters 

Normally the boundary of a particle trajectory plot 
is adjusted to correspond to the boundary of the highest grid 
into which any of the N emitted particles fell during tracking. 
Unfortunately particles circling in small magnetic fields 
tend to make very large loops, i.e., into grid 12 or higher. 
Thus the automatic plot boundary selection can obscure de- 
tails of individual particles returning to the object in 
grid 1. (The object silhouette is only plotted if the 
boundary of the plot corresponds to grid 6 or lower and fine 
details are usually only apparent if the boundary of the plot 
is set to grid 4 or lower.) In addition there may be some 


cases in which it is desirable to have the same set of 
particle trajectories plotted more than once, each time 
using a different maximum grid boundary for the plot* The 
following three keywords provide the user with the capability 
to deal with these problems and special requirements, 

NGBND Grid number to which the boundary of the particle 
trajectory plots corresponds to. (Note that all 
grids which are higher than 2 but fall inside the 
plot edge are automatically outlined for reference 
on the trajectory plots.) Acceptable range is 
integral NGBND > 0 , i . e . , NGBND « 0 , 1 , 2 , . . . 

If NGBND » 0 then the boundary of the plot will be 
automatically adjusted to correspond to the highest 
grid into which any of the N particles fell during 
tracking. Default value is NGBND "0. 

NGPLOT Number of different plots of same set of trajectories 
desired (only has an effect if NGBND ^ 0) # Each 
trajectory plot is made using a different grid bound- 
ary (see NGINC) . The acceptable range is integral 
NGPLOT such that 1 £ NGPLOT £ 4 . The default is 
NGPLOT = 1. 

NGINC Increment factor for grid boundaries of successive 
trajectory plots. If NGPLOT > 1 then the same tra- 
jectories will be drawn on NGPLOT different plots. 

The grid boundaries are incremented for each succes- 
sive plot as follows: 

IG = NGBND + (I-l) t NGINC for I * I, . . . , NGPLOT 
Acceptable range is integral NGINC £ 0. Default 
is NGINC =0. 

Note that the highest grid value selected for plotting 
boundaries by NGBND, NGPLOT, and NGINC should not exceed LIMGRD 
if useful plots are to be obtained. 
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MOOPAR 


LIMGRD 


ZETA 


ELMASS 

PRM.ASS 

ELCH 


Modulo number for particle trajectory line segment 
plotting (significant only if PLPART option active) . 
MODPAR > 0 results in every MODPAR*^^ trajectory 
point being plotted, MODPAR « 0 or keyword option 
omitted will result in the default of every 2 ^^“^ 
point being plotted where IR is the index number 
of the NASCAP computational grid in which the 
particle is located. Thus the default is for 
dynamic adjustment of MODPAR. 


Highest grid in which particle tracking is permitted. 
If particle passes outside this grid it will be 
assumed to have escaped to infinity. Default value 
is LIMGRD »< NG, the highest NASCAP computational 
grid. Acceptable range is LIMGRD ^ 1. (Note that 
particles which exit from the highest NASCAP com- 
putational grid NG are tracked using a monopole 
potential.) 

Angle of particle detector on a boom cell. (Default 
is ZETA = 0.0“.) Not relevant if ICELL is the 
index of a surface cell of the satellite. 


/ If th6- user intends to modify particle mass and/or 
I charge using the following four keywords then he 
must be certain that the same particle masses and 
charges are specified on the environment FLUX 
([^definition file! 

Electron mass used for reverse trajectory tracking. 
(Default is ELCH » 9.1091 x lo“^^ kG.) 




Proton mass used for reverse trajectory tracking. 
(Default is PRGH = 1.67252 x lo"^^ kG.) 

Electron charge. (Default is ELCH = -1.602 x 10 
coul. ) 




mrnsmmmm. 


PRGH Proton charge. (Default is PRCH ■ +1.602 x lo”^^ 

coul.) 


Description of Old Detector Keyword Options Which Have Been 
Mo<ii fied : 


ICBLL 

NSTP 


DE 

DEK 


VCODE 


N 

PSCALE 


Now may be index number of a boom surface cell. 

Maximum numt<er of steps per particle allowed during 
reverse trajectory tracking. Acceptable range is 
1 < NSTP < 30,000. Default is NSTP - 500. As a 
rule of thumb, assuming VCODE ■ 0.3, NSTP should 
be increased by at least 400 for each additional 
NASCAP grid used. If NG - 3 then NSTP ■ 1200 is 
probably not unreasonable. 

Deleted keyword. Function superseded by DEK. 

Width of detector energy aperture in percent. 
Acceptable range is DEK > 0.0. Default is DEK > 0.0. 
If the energy the detector is looking at is E then 
the flux value will be obtained from Eq. (8.9), 
p. 156, of Reference 4, where 

( DEK \ 

1 (constant % energy window) . 

If ICELL is the index number of a boom surface cell 
then the value specified for VCODE is automatically 
scaled by where IG is the index number of 

the grid in which the boom cell is located. 

Definition of this keyword is same as before but 
default value has been changed from 100 to 20. 

Definition of this parameter is same as before 
except that documentation incorrectly stated that 
the proton flux multiplication factor was 10**PSCALE 
when in fact it is just PSCALE. The default is 
PSCALE = 10^. 


2.6 EMITTER ROUTINE KINENG 
Definitions ; 

- eleotro.tstio potenti.l o£ emitter conductor 

■ electrostatic potential at top of dielectric layer covering 
conductor 

Eq ■ particle kinetic energy at conductor surface 

6^ » polar emission angle of particle at conductor surface 

■ particle kinetic energy at top of dielectric layer covering 
conductor 

■* 

6, ■ polar emission angle of particle at top of dielectric layer 
covering conductor 

m « particle mass 

q * particle charge 

V , v„ particle velocity components at conductor surface 
X y » 

v^f v^ particle velocity components at top of dielectric 
X y 

layer covering conductor 
Motivation ; 

Particles are conceptually emitted from the geometric 
center of the conductor underlying the surface cell on which 
the low-density NASCAP particle emitter is located. In NASCAP 
potential nodes which lie on surface cells are considered to 
represent potentials on the surface of the thin dielectric 
coating which overlays the surface cell conductor rather than 
the potentials of the conductor itself (see Figure 2.2). 
Therefore, it is not possible to numerically solve the equa- 
tions of motion for a particle in the dielectric coating region 
using electric fields obtained from the NASCAP routine EFIELD. 
(EFIELD does not know about the dielectric region.) NASCAP 
emission routine EMITER circumvents this difficulty by 
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employing an approximate analytic treatment within the di- 
electric region. 



Figure 2.2. 2-D illustration of surface cell/dielectric 

interface. Shaded area is dielectric region. 


Analysis ; 

The analytic treatment of particle motion in the di- 
electric layer covering an emitter surface cell conductor 
is made using the following set of assumptions: 

1. Particles in the dielectric region are accelerated 
only by an electric field which is normal to the plane 
of the conductor . 

2. The point on the dielectric surface where the particle 

passes into free space lies on that line which is 
normal to the plane of the conductor and passes through 
the geometric center of the surface cell. (This is 
an approximation.) . 
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Energy conservation applied to Figure 2.3 yields the 


relation between initial particle kinetic energy E and 


emission polar angle 6^ at the conductor surface and the 


kinetic energy E^^ and polar emission angle of the particle 
at the top of the dielectric layer as follows: 


Eq - m/2(v^ + v^) 


Ej^ * m/2 


= m/2(v'2 w v/) 
or rewriting 


(initial kinetic energy) 
(final kinetic energy) 


(since ~ ^ consequence of 

normal electric field) 


2 2 
X m y 


But 




2 

sin^ 9 
in o 


so 


V 


.2 ^ 2 I 

X ~ m \ 


. 2 

Et - sin 0^ 
L O O 


( 2 . 1 ) 


Also 


m V. 


.2 


cos 0- = 


L 2 E, 


( 2 . 2 ) 


and 


Et = E^ + AE 
L O 


(2.3) 


where AE = q(c!)_ - d>_ ) 

C li 
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Center 




Conductor 


Dielectric Layer/ 
Vacuum Interface 


0 1 

0 < < 7t/2 


Figure 2.3. Initial and final particle velocity for case 
AE > 0. 


Putting Eq. (2.3) into Eq. (2.1) gives 


v^2 = 2 ( 

X m \ 


- - ( E 
m \ o 


E^ + AE - E^ sin 0^ 
o o o 


cos 6^ + AE 
o 


Then using this in Eq. (2.2) we obtaini 


, E^ cos 0^ + AE 

^ ^ o o 

®L ' E^ + AE 


0^ = cos 




cos 9^ + AE 
o o 








Implementation ; 

The analytic treatment presented above is implemented 
in the NASCAP particle emitter routines via subroutine 
KINENG. Given the emission energy and angle 9^ at the 
conductor surface, KINENG returns the particle energy E^ 
and angle 0^ at the dielectric/vacuum interface as given by 
Eqs. (2.3) and (2.4) . 

In the event that Eq. (2.3) yields a negative energy or 
Eq. (2.4) has no real solution, KINENG sets a flag to indicate 
that the particle can never reach the dielectric surface and 
hence its emission is suppressed. 

Note that an acceptable solution of either equation 
alone is not sufficient to guarantee that a valid solution of 
the other equation exists! 

Examples: 

1. E- = E /2 satisfies Eq. (2.3) but (2.4) has solution 

L O 

only for 0^ < 45® . 

2. 0^ = 0 •+ = 0 for any AE, but some AE will make 

Ej^ < 0 violating conditions on Eq. (2.3) . 


Operational Note ; (add to p. 182 of Reference 4) 

If the NASCAP user includes the IPRNT option card in 
his particle emitter keyword file then prior to the one-line 
summary for each particle trajectory completion KINENG will 
print one of two possible informative messages pertaining 
to the trajectory. The two possible messages are of the 
form 

1. ***EMISSION SUPRESSED BY KINENG*** 

is printed if Eq. (2.3) and/or (2.4) had no valid solution. 

2. **IN KINENG THETA = CHANGED TO 

VANOD = CHANGED TO — 


Message 1 may appear up to a maximum of 50 times. If this 
limit is exceeded a warning is issued and no further such 
messages will be printed. 

The variables in these messages coXi.espond to the 
definitions as follows: 

THETA = 8^ (radians) PAVG * <1)^ (volts) 
o u 

VANOD = E (eV) VC = <j> (volts) 

o c 

The "changed to" variables are 9^^ and Ej^. 

If message 1 is not printed by KINENG it is still 
possible to receive another similar message: 

♦♦EMISSION SUPHESSED BY LOCAL PIELD^^ . 

This message is generated by the EMITER routine in response 
to a flag which is set within the PUSHER routine. The 
PUSHER routine contains a preliminary test to determine 
(based on the emitter local surface cell field) whether or 
not the particle has sufficient energy to travel a distance 
of more than one cell. If it doesn't then there is no point 
in explicitly traclcing such a particle. 

2.7 BOOM CELL MODELING 

(This section describes a feature that has been in 
NASCAP for some time, but it was not included in last year's 
Final Report . ) 

To treat the boom one constructs a mixed geometry 
(cylindrical-cartesian) finite element solving technique, 
starts by enclosing each section of the cylindrical boom in 
four cubes of sides of unit length, which are part of the 
cubical 3-D mesh. It is apparent that it is enough to treat 
only one quadrant. 
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The object is to construct the finite element matrix 
that represents the finite element interpolation potential 
in the region of Figure 2.4. 

The interpolation potential has to have the following 
properties: 

a. It has to be constant on the circular section of 
the boom. 

b. It has to satisfy Laplace's equations in the 
region of Figure 2.4. 

One such interpolation potential is 

V(r,(J)Z) « A(r,(j)) (1-z) f B(r#(^)z , 

where the coordinates system is that defined in Figure 2.4 
with 

A(rf(j)) ~ 01^ "t" An ^ ^2 

and B(r,(J)) is the same functional expression with a's replaced 
by 3's. • 

We call 32^ the potentials on the three avail- 
able corners of the square (counting counterclockwise) and 
a^ the potential over the quarter circle (see Figure 2.4) on 
the plane z = 0 , and on the other hand bj^ , b 2 / b^ the equiva- 
lent entities on the plane z = 1. Then the a's bear a linear 
relationship with the a's obtained by asking A(r,cj)) to take 
the values a^, a^^, 32# a^ on the quarter circumference and 
the three corners respectively and equivalently for the $'s 
and b's. 




Then we construct the 8x8 matrix, to be used as a 
submatrix of the total mesh, such that 

f |VV|^ rdr d(f> dz - X S X , 

where 

•» 

X “ (a^ • • • f t • • • f ^4^ 
and X is the associated column vector* is given by 

' 0 1 i J 
' 7 1 ^ 1 J 

We have written a subroutine to obtain S as a function of R 
with arbitrary accuracy, because we have done the quadratures 
needed analytically. 

Figures 2. 5-2. 8 present potential contours resulting 
from a two-dimensional calculation on a 9 x 9 mesh contain- 
ing a central boom and a grounded outer conductor having 
diameter equal to the mesh diagonal. The boom radii range 
from 0.01 to 0.32 mesh units. In the figures the boom out- 
line is represented by the second innermost contour. With 
increasing boom size the fields at the boom surface decrease, 
and the fields in the surrounding space increase. The cal- 
culated boom capacitances are one to two percent larger than 
the theoretical values. Distortion of equipotentials near 
the boom cell boundaries results from the differences in 
interpolation functions between boom and bilinear cells. 




X i sin? 


0 < z < 1 
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3 . USABILITY 


NASCAP is constantly being improved to fit the needs 
of users. No matter how accurate a computer program may be, 
it is only useful to those who can run it. The changes we 
have made to make NASCAP more usable, fall into three cate- 
gories: input, output, and new features. 

"Input” includes all those changes for which NASCAP 
still does the same thing, but with far less persuasion and 
patience on the user's part than was previously necessary. 
First, these changes reduce problem set-up time. They allow 
users to concentrate on physics more, and keypunching less. 
Second, they save computer time by reducing the number of 
wasted problem runs necessary to produce a good one. 

"Output” changes improve the selection of information 
available to the user. Most of the changes take information 
that NASCAP already stores internally, and make it available 
as printout or graphics plots. One of the changes (print 
flags) also concerns some information that NASCAP already 
prints out. It lets the user say, "Keep it to yourself, 
NASCAP. I don't want to hear about it." 

The "New Features" in this section (Deadline and Timer 
Spin Command, and Automatic Convergence) are some of the most 
widely asked-for improvements. In their various ways, they 
also save computer time and simplify the job of running 
NASCAP . 
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3.1 


INPUT 


While the NASCAP code is constantly becoming more 
sophisticated, it is at the same time becoming easier to use. 

As the number of users grows, the new users suggest ways to 
simplify input procedures. 

Perhaps the biggest procedural change made in the last 
year is the new free format input. This allows users to ignore 
column numbers and format statements. 

Another big change is the introduction of default 
material properties. NASCAP now provides a set of material 

I 

properties automatically. 

The user input for defining booms and the input for 
special cell selection have been improved. 

In all of these cases, the old input procedures work 
as before. The new, simplified procedures offer an addi» 
tional alternative. 

Finally, there is a new way to specify the flux pattern 
for a test tank electron beam. The old method (which still 
works) requires a polar coordinate pattern. The new method 
uses a beam pattern specified in cartesian coordinates. 

Appendix A of this report is a replacement of Section 7.3 
(Flux Definition File) of the NASCAP User's Manual. It de- 
scribes both techniques in full. 

3.1.1 Free Format Input 

This year's NASCAP has incorporated a Pree Format 
feature. This means that you can scatter your input lines 
right across the page, without paying any attention to column 
numbers. 

No longer will NASCAP change your number of computa- 
tional grids from 2 to 20, just because you put your 2 in 
the wrong column. You can now ignore the Format statements 
in the NASCAP User's Manual. 


There must be a space or a comma separating each input 
value, but otherwise NASCAP doesn't care how close together 
or far apart they are. The card images: 

SUNDIR 2,2,1 

SUNDIR 2 21 

are equivalent. 

There are still a few input restrictions. You must 
put the right number of items on each card. You can't leave 
blanks to indicate zeroes. And don't put a decimal point on 
a number that NASCAP expects to be an integer. The forms 
"2" and "2." are both acceptable for REAL number (floating 
point) input. But for INTEGER (fixed) numbers, only "2" 
would be acceptable. 

If NASCAP doesn't like your input, then subroutine 
FREED or subroutine CLASFY will probably print out angry 
error messages. 


3 . 1 . 2 Surface Cell Specifications 

NASCAP prints a flux breakdown and other information 
for surface cells that are specified by the user. There are 
two ways to specify which surface cells are of interest, the 
SURFACE CELL card and the SURFACE AT card. Both cards go in 
the "options" file to be read at RDOPT time. 

The SURFACE fSLL card has one integer — a surface cell 
number. This number can be obtained by finding the desired 
cell on the surface cell list printed by OBJDEF. 


The SURFACE AT card has three or six integers to de- 
scribe the geometrical location of a surface cell. The 
format is 


or 


SURFACE 

AT 

ix 

jy 

kz 

SURFACE 

AT 

ix 

jy 

kz 


ny 


nz , 
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where (ix,jy,kz) is lower-left corner (in absolute coordinates) 

of the volume cell in which the surface is located or out of 

which it points, and (nx,ny,nz) (optional) is the surface 

normal. Some examples are shown in Figure 3.1. Note than 

an invalid cell specification is interpreted as cell 1251, 

so as not to interfere with validly specified cells. 

rupfaoe at -iO -10 IS 
SURFACE AT 10 10 IS 

FNDCEL FOUND CELL NO. 238: 111212176041. 

CELL 288 INFORMATION TO BE PRINTED. 

surface at 10 10 16 
SURFACE AT 10 10 16 

FNDCEL FOUND CELL NO. 289: 011212206001. 

CELL 239 INFORMATION TO BE PRINTED. 

surface at 9 9 14 
SURFACE AT 9 9 14 

FNDCEL FOUND CELL NO, 244: 011111160101. 

CELL 244 INFORMATION TO BE PRINTED. 

surface at 4 9 17 
SURFACE AT 4 9 17 

FNDCEL 000411210000 NOT FOUND 

CELL 12S1 INFORMATION TO BE PRINTED. 

surf ace cell 48 

QlIPPiiPF TFI I AQ 

CELL 43 INFORMATION TO BE PRINTED. 

surface cell 49 

f^lIRPArP' PFI I 49 

CELL 49 INFORMATION TO BE PRINTED. 

surf ace ce I I 13 
SURFACE CELL 13 

CELL 13 INFORMATION TO BE PRINTED. 

Figure 3.1. Examples of "Surface Cell" and "Surface At" keyword 
options. An invalid cell specification is inter- 
preted as cell 1251. 


3.1.3 Boom Input 

A new form of input is available for the AXIS card in 
BOOM definition in NASCAP, The new form is (free format) 
AXIS IBX IBY IBZ lEX lEY lEZ 
where (IBX, IBY, IBZ) and (IEX,IEY,IEZ) are the beginning and 
ending locations of the boom axis, respectively. The loca- 
tions are all assumed to be in first grid coordinates using 


this form (even if the boom extends outside the first grid) . 
Any offset is applied to both beginning and ending locations. 
Locations in outer grids must^ of course, correspond to nodal 
points. The old form of the AXIS card is still valid; 

AXIS IBX IBY IBZ IBG lEX lEY lEZ lEG 
where IBG and lEG are the grid number of the starting and 
ending axis locations, respectively, and where any offset is 
applied only to locations in the first grid. 

An example of the new input form: a boom extends from 

the center of the inner mesh along the negative X axis just 
one node into the third grid. The default offset grid center 
is used. 

New form: 

AXIS 0,0,0 -20,0,0 

Old form: 

AXIS 0,0,0 1 4,0,0 3 

Both axis cards will describe the same boom. NASCAP 
knows which form you are using by the number of integers on 
the AXIS card. 

3.1.4 Default Material Properties 

Default values for material properties of twelve mate- 
rials have been provided in NASCAP. A list of the default 
materials and parameters is given in Appendix B. The routine 
MATDEF contains the default values in the code. If proper- 
ties are not read in for a material during object definition, 
the first reference to that material name will cause the 
default values to be used. The default values will be 
ignored whenever a set of properties is specifically read in 
for a material in the default list. 

A new command has been added to NASCAP to allow mate- 
rial properties to be changed without redoing OBJDEF. The 
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conunand card 'NEWMAT LUN' ^ where LUN is an integer, causes 
revised material parameters to be read from unit LUN. For 
each material, LUN should contain a card with the material 
name, followed by three cards with the revised parameters 
exactly as during OBJDEF. An 'END* card terminates reading 
from LUN. The NEWMAT segment reprocesses material para- 
meters, recalculates the large caivacitances between sur- 
face cells and the underlying conductors, and recalculates 
the conductivity matrix elements. There is one restriction 
on parameter changes: the sign of parameter 3 (bulk conduc- 

tivity) and 14 (surface conductivity) cannot be changed. 
Therefore, a material initially treated as a perfect con- 
ductor (parameter 3 negative) cannot subsequently be treated 
as an insulator (parameter 3 positive) , or vice-versa. Also, 
a material for which surface resistivity was ignored (para- 
meter 14 negative) cannot subsequently be assigned a surface 
resistivity (parameter 14 non-negative) , or vice-versa. 

NEWMAT may be called any time after OBJDEF. 

3 . 2 OUTPUT 

The major change in NAS GAP output has been the develop- 
ment of post processor codes through which the user can inter- 
actively request specific information. However, the basic 
NASCAP output has been improved in four areas — current balance 
printout, flux breakdown, print flags, and contour plot flexi- 
bility. 

The current balance printout shows, for each conductor, 
where the current is coming from. The flux breakdown does the 
same for surface cells, with some additional information about 
cell geometry and fields. The new print flags allow a user to 
control the amount of NASCAP printout. And the contour plot 
flexibility feature lets the user tell NASCAP how many grids 
should be shown in potential contour plots. 


3.2.1 Current Balance Printout 


An analysis of current balance to conductors is pro- 
vided at the end of the LIMCEL section. For each conductor, 
the quantities provided are; 

a. Net current (implicit) 

b. Plasma current (explicit, taking account of 
limiting of low-energy electrons) 

c. Conductivity current (implicit) 

d. "Battery” current = (Net current] - [Plasma current] 

- (Conductivity current] 

In steady state, the "battery" current represents current 
supplied to maintain biased or fixed potentials on conductors. 


3.2.2 Flux Breakdown 


The flux breakdown printout has been modified to dis- 
play more information, as shown below: 


SURFACE CELL NO. 237 


CODE = 011310050302 

DOCATICN = 11 8 5 

NORMAL = 00-1 

MATERIAL = SOLAR 


POTENTIAL = 7.441+00 VOLTS 
STRESS = -1.067+04 VOLTS/METER 


EXTEmL FIELD = 

FLUXES IN A/M**2 

INCIDENT EILECTRCNS 

RESULTING BACKSCATEER 
RESULTING SECCNDARIES 
INCIDENT PRCTCNS 

RESULTING SECCNDARIES 
BUU^ CCNDUCTIVITy 
PHOTOCURRENT 

NET FLUX 


8.867+00 VOLTS/METER 

1.98-07 

8.48-08 

1.30-07 

7.25-09 

5.29-09 

-1.07-13 

6.04-06 

6.06-06 


Figure 3.2. Flux breakdown printout. 
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STRESS is the internal field for insulating surfaces 
(previously labeled 'FIELD' in the printout) . The external 
field was not printed before; this value is used (when posi- 
tive) in the field limiting of low energy emission. The 
initial bulk conductivity is now also printed for insulating 
surface cells. 

To request a flux breakdown printout for any surface 
cell, the user inserts a card in the options file, read by 
RDOPT. There are two forms for this card — one is 'SURFACE 
CELL n' where 'n' is the surface cell number. The other is 
'.SURFACE AT i,j,k' where 'i,j,k' is the cell location. Both 
are described in Section 3.1.2, "Surface Cell Specifications". 

3.2.3 Print Flags 

Keyword input may be used to vary the print level 
from various portions of NASCAP, 

'PRINT POTENT' results in an increased level of out- 
put from the large Scaled Conjugate Gradient potential solver. 
This is intended mainly for debugging purposes. 'NOPRINT 
POTENT' restores the default. 

'PRINT OBJDEF' results in printing of the element 
table and some point-lists which are now by default suppressed, 
'NOPRINT OBJDEF' suppresses printing of the surface cell list 
and list of material properties. 

'PRINT LIMCEL' causes the results of the preliminary 
charging analysis to be printed prior to entering the final 
charging analysis. For each insulating surface cell and con- 
ductor the initial and trial potentials, external field, cur- 
rent level, and status (fixed or not fixed) are given. 

'NOPRINT LIMCEL' causes reversion to the default. Also, 
some obsolete debugging prints have been removed from the 
LIMCEL section. 







3 . 2 . 4 Contour Plot Flexibility 

A new keyword, NGPLOT, has been added to NASCAP's 
RDOPT routine. This allows the user to specify the number 
of grids to be considered in potential contour plots. Using 
the keywords ICON, NCON, and NGPLOT, contour plots will be 
generated as indicated below: 


ICON 

NCON 

NGPLOT 


0 

0 

any 

None 

0 

any 

0 

None 

>0 

0 

0 

Three views each of 1 and NG 
grids when ICYC = 1 mod ICON. 

0 

>0 

>0 

Extra views only plotted in 
NGPLOT grids every timestep. 

>0 

>0 

>0 

Three default views plus extra 
views plotted in NGPLOT grids 
every timestep. 

>0 

>0 

0 

Three default views plus extra 
views plotted in 1 and NG grids 
when ICYC = 1 mod ICON. 

>0 

0 

>0 

Three default views plotted in 
NGPLOT grids when ICYC = 1 mod 
ICON. 


Figure 3.3. Keyword values for contour plots . 

Note that additional flexibility is available through 
use of the interactive routine POTPLT. 
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3.3 


NEW FEATURES 


Three new features in NASCAP for the 1980 's are; dead- 
line and timer, automatic convergence monitor, and the spin 
command . 

The deadline is quite useful in allowing for scheduled 
computer down time. Rather than worry every night that the 
computer will dump you before your NASCAP run finishes, you 
can set a NASCAP deadline, after which the program gracefully 
terminates. Timer shows the user how much execution time is 
used by various parts of NASCAP. 

The automatic convergence monitor eliminates another un- 
certainty for the NASCAP user. It eliminates the old keyword 
MAXITR, which was hard to figure out, and was also critical to 
the successful running of NASCAP. 

The SPIN command will save much computer time analyzing 
satellites that turn in relation to the sun. It is often de- 
sirable to run NASCAP timesteps of ten minutes or more. But 
a spinning satellite used to require several timesteps per 
rotation to properly figure solar illumination. Now, an 
average illumination can be calculated. 

The SPIN command is one example of a code change that 
appeared to be quite simple, but ended up causing changes in 
many parts of the code, including the hidden line routines 
(HIDCEL) and the longtimestep routines (LIMCEL) . 

3.3.1 DEADLINE AND TIMER 

The keyword specification 

DEADLINE hhmmss 

will cause TRILIN to exit rather than start another timestep 
after the host computer's internal system clock has passed 
through the DEADLINE time. 

Examples : 

DEADLINE 235959 


will not allow a timestep to start after midnight of the day 
the deadline card is read. 

DEADLINE 060000 

will not allow a timestep to start after 6:00 a.m. (the next 
day, if the DEADLINE card is read during the evening) . 

Remember that the effect of the DEADLINE is gone after 
exiting (and re-executing) NASCAP. 

The TIMER option gives the user a better idea of how 
NASCAP execution time is used up. If you insert a card 
reading "TIMER” in the options (RDOPT) file, then NASCAP 
will print out the accumulated amount of execution time, at 
various points in the calculation. 

3.3.2 Automatic Convergence Monitor — MAXITR 

It is no longer necessary for the NASCAP user to guess 
how large to set MAXITR to achieve convergence o:: the potential 
solver. MAXITR can be totally left out of the options file. 

NASCAP now keeps track of the measure of convergence and 
its rate of decrease. The iteration procedure terminates if 
the following three conditions are satisfied; 

1. The number of iterations is at least 8* (no. of grids) . 

2. The measure of convergence is decreasing. 

3. The measure of convergence has been reduced POTCON 
orders of magnitude from its initial value. 

The parameter POTCON is set by default to 8 in CAPACI and 
to 4 in TRILIN. It may be changed by keyword input in the op- 
tions file. As a "measure of convergence" in the potential iter- 
ations, the square of the length of the residual vector is used 
(i.e., p" in Figure 7.1 of the NASCAP User's Manual^^^). Since 
MAXITR will' now not be normally used to terminate the convergence 
process, its default value has been changed to 99. 
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3.3.3 SPIN Command 

A SPIN coiranand has been added to the main NASCAP routine. 
This command generates a table of averaged solar illumination 
for each surface cell of a rapidly spinning satellite. SPIN 
performs a straightforward numerical average of solar illumina- 
tion factors for a sequence of sun directions about a specified 
rotation axis . 

The SPIN routines require specification of NVIEWS, 
the number of sun directions to include in the averaging ^ 
and SPINAX(3), a vector specifying the rotation axis. The 
default values are NVIEWS = 8 and SPINAX = (0.0.1.); these 
values are used when LUN * 0 is supplied on the SPIN com- 
mand card. If LUN = 5, SPIN reads two cards from the input 
runstream. For any other positive value the cards are read 
from the specified file. The first card read contains NVIEWS 
(an integer) and the second card contains the x, y, and z com- 
ponents of the spin axis. The spin axis vector need not be 
normalized. Tests with NVIEWS = 8 indicate that the numerical 
averaging procedure is accurate to within a few percent for 
simple objects. 

The SPIN driver routine generates NVIEWS calls to 
HIDCEL/ with the sun initially at the original input sun 
direction/ and then rotated uniformly thereafter. If the 
CONVEX option is specified, calls to HIDCEL are not required. 
Once the SPIN command is invoked, NASCAP assumes all later 
TRILIN Cycles (including restarts) for the satellite will use 
the resulting tables of averages illumination, unless HIDCEL 
is explicitly called at a later point. 

Figures 3.4 and 3.5 illustrate the equilibrium poten- 
tials reached using the SPIN option for a one-grid SCATHA 
model . 
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Figure 3.4. Potential contours around satellite at equili- 
brium in sunlight;# with conductor potential at 
-700 volts. Average solar illumination applied 
with initial sun direction = (.1,1. /I.) and spin 
axis = (1. ,0. ,0.) . 








4 . CODE REVIEW 


In a large sense, all of the changes made to NASCAP 
are a result of code review. This section of the report is 
devoted to those changes which are general or miscellaneous 
or don't fit under "Physics”, "Usability", and "Auxiliary 
Codes". The six categories here are: centering, faster 

interface routines, 12150 surface cells, code cleanup and 
commenting, shadowing (HIDCEL) improvements, and 'TYPE* con- 
flict. 

The new NASCAP centering scheme was a major improve- 
ment. It increased the stability and accuracy of the code, 
while at the same time greatly reducing running time! Faster 
execution time was also achieved by speeding up the subrou- 
tines that do calculations on the NASCAP grid interface 
regions . 

The number of allowed surface cells was increased fzoni 
1024 to 1250. A comprehensive code cleanup made the Fortran 
coding easier to read and understand. Some improvements in 
the shadowing routines (HIDCEL) removed limitations on some 
unusual shadowing situations. Finally, we removed a small 
but annoying problem that resulted from conflicting flux type 
specifications. 

4.1 CENTERING 

The centering scheme used to update surface charges 
and potentials has been completely rewritten, a major task. 
Cell centered potentials and charges are now stored for each 
insulating cell and for each conductor. (These cell poten- 
tials are printed every timestep.) Potentials at nodal 
points bordering disparate materials are now fixed to area 
weighted averages. This scheme, exactly the reverse of the 
previous centering, results in much smoother potential 
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variations through transition regions. This scheme also 
allows the use of the capacitor model to fix surface poten- 
tials while space fields are converging, which significantly 
improves the rate and stability of the iterative convergence 
in POTENT. 


0 

0 

' 1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Figure 4.1. Centering scheme. Dark spots are nodal points — 
open circles are cell centers. Cell potentials 
and charges are now calculated at cell centers 
rather than at nodal points. 


4.2 PASTER INTERFACE ROUTINES 

The interface routines called by COPROD have been re- 
vised to make use of fewer subroutine calls and to decrease 
the number of matrix operations required. As a result of 
this improvement and the recentering scheme, a full run of 
the four-grid SCATHA model to equilibrium in eclipse now 
requires approximately two hours UNIVAC 1100/81 CPU time, 
whereas nearly eight hou„i« were required priviously. 



4.3 


1250 SURFACE CELLS 


NASCAP is a large code. It always nudges the size 
limitations of its host computer. One result of this size 
problem is that NASCAP objects have always been limited to 
1024 total surface cells. 

This year, by using sophisticated mapping techniques 
and shaving space wherever possible, we have increased to 
1250 the number of surface cells allowed. Because of extra 
arrays needed for insulating surface cells, not more than 
1024 of the 1250 total are allowed to be insulators. 

4.4 CODE CLEANUP AND COMMENTS 

Dr. M. J. Mandell went through the entire NASCAP code, 
subroutine by subroutine. He performed a major code cleanup 
and added comments. In detail; 

1. The SUBROUTINE or FUNCTION statement was placed 
at the head of each symbolic element. 

2. A comment was added stating the calling routine 
for each subprogram and the code module to which 
the subprogram belongs . 

3. When appropriate, references were added to NASA 
CR's 159417, 135259, 159595, and 135256. In a 
few places references to open literature were 
added. 

4. Excess lines of coding were deleted. Such lines 
included multiple spacer (blank comment) lines, 
debugging print statements , unused FORMAT state- 
ments, unnecessary COMMON statements, and old 
algorithms which have been superseded. 

5. In a few places loops were replaced by S3MOVE 
or S3ZERO. 
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6. Obsolete subprograms were expunged. 

7. Diverse sorting routines were replaced by a 
single, faster subroutine, SRTWDS. 

4.5 SHADOWING (HIDCEL) IMPROVEMENTS 

Several improvements have been made in the operation 
of the HIDCEL shadowing routines. With the operation of the 
ROTATE and SPIN features of NASCAP in conjunction with the 
use of long thin booms, some difficulties were encountered 
in the shadowing algorithms, especially in the treatment of 
small boom "tips” in the SCATHA models. These operational 
problems have now all been corrected. Figure 4.2 illustrates 
one such shadowing case which caused the code to fail: it 

contains A2 (shadow ing ) surfaces which are entirely contained 
within A1 (shadowed) surfaces. The HIDCEL algorithm now 
recognizes the occurrence of such cases and calculates the 
correct shadowing fraction by simply subtracting the A2 area 
from the A1 area. Another improvement involved a new routine, 
POLYN, which eliminates redundant vertices from generated A1 
polygons as necessary; previously HIDCEL could malfunction 
when polygons having over twelve vertices were generated. 

4.6 'TYPE' CONFLICT 

In the original NASCAP, a conflict between the flux 
type in the flux file and the type chosen in the options 
(RDOPT) file would halt program execution. The user was 
required to supply information that was either redundant 
or incorrect. 

Now, the flux type specified in the flux file is the 
one that counts. The options (RDOPT) file does not need to 
specify what flux type is being used. 



5 . AUXILIARY CODES 


NASCAP was designed as a stand-alone code. It re- 
quires no other programs for problem setup or interpreta- 
tion. But now there is a small family of associated codes 
that make NASCAP use easier/ and give the user a new degree 
of flexibility. The major auxiliary codes are TERMTALK/ 
MATCHG/ and the workshop codes (FILES/ OBJCHECK/ and POTPLT) . 

In original NASCAP/ all output comes back as hard 
copy printout or graphics. If you want to see the evolution 
of electric charging on a particular surface cell/ you have 
to page through a stack of output and read the potential at 
each timestep from a list. TERMTALK gets the computer to 
sort through NASCAP restart files to find the data you 
want. 

MATCHG is a special code for investigating material 
properties. In effect, you take a little patch of any sur- 
face material, and subject it to a specified flux. You find 
out how that material will charge without the influences of 
geometry/ other materials, and sunlight. This gives you a 
general idea of what will happen to a test object. MATCHG 
is very much faster and cheaper to run than NASCAP. 

The workshop codes were designed to simplify things 
for new NASCAP users. They were used for the NASCAP work- 
shops at NASA/LeRC. FILES does Univac file manipulation. 
OBJCHECK draws object illustration plots. POTPLT draws 
electric potential contours. 






5 . 1 TERMTALK 

TERMTALK is an interactive program for retrieving 
NASCAP data. It accesses NASCAP restart files to produce 
charging history tables and graphs and a variety of other 
information. 

TERMTALK is a menu-style program - the user requests 
the desired information rather than answering a series of 
yes or no questions, "Do you want such and such? How about 
such and such? ...” As a result, the user can get any of 
the types of available output just by typing one or two com- 
mands, no matter what section of the program is currently 
operating. 

There are four main program modules, named HISTORY, 
LATEST, SINGLE, and SPECIAL. 

The HISTORY module outputs for any surface cell the 
time history of the user's choice of five quantities — 
potential, electric flux, external electric field, internal 
electric field (stress) , or potential difference between an 
insulator surface and underlying conductor (delta) . The out- 
put comes in the form of a printed table, and/or a rough plot 
of quantity versus time. The printed table will show up to 
seven surface cells across the page. 

The LATEST module gives the user a complete list for 
all of the surface cells and conductors of any of the five 
quantities mentioned above. They may be printed in sequence 
or ordered by magnitude.' Partial listings, including only 
some surface cells, are available. 

Module SINGLE prints out information about a single 
surface cell. In addition to the five standard quantities, 
you can get the cell location, its surface material, its 
shape, its normal, and the number of its underlying conductor 
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The SPECIAL module lets the user "turn off" all output 
coming to the terminal. This unseen output can be printed 
on a line printer at the end of the session. SPECIAL also 
allows the user to change the NASCAP cycle number. 

PROGRAM STRUCTURE 

TERMTALK is organized into four modules — LATEST, 
HISTORY, SINGLE, and SPECIAL. Each of the modules has its 
own set of commands. Each module also has access to a set 
of "aid" routines, named HELP, SUBSET, OUTLINE, LOCATION, 
and EXIT. A crude diagram of the program is shown in Fig- 
ure 5.1. 

The user moves from any module to any other by typing 
the name of the new module. "Aid" routines always return 
to the module from which they were called. 

CELL SPECIFICATIONS 

NASCAP assigns a cell number to each surface cell. 

Each cell can also be identified with five basic geometrical 
attributes — location, material, surface normal, shape, and 
conductor number. TERMTALK provides a way for the user to 
go back and forth between cell numbers and geometrical at- 
tributes. 

'SINGLE' takes a cell number and prints out attributes. 
It also prints out charging information. 

'SUBSET' takes the user's geometrical specifications, 
and determines the set of cells that satisfy these. It is 
extremely general. You can easily determine the number of a 
particular cell, or you can define a complicated group of 
cells. For example, you could define the group of all non- 
boom cells that are either KAPTON or TEFLON and lie on con- 
ductor #1; or all cells between the planes of X = 1 and 
X = 5. 
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Figure 5.1. Diagram of TERMTALK. 







CHARGING HISTORY 

TERMTALK can plot surface cell potential as a function 
of time. It also prints cell histories in a tabular form. 

And it gives tables and plots of other quantities such as 
electric flux and electric fields. 

If you are more interested in final equilibrium values 
than in histories, the module ’LATEST' gives only the most 
recent information. 

INSTRUCTIONS FOR USE 

‘ The most extensive documentation for TERMTALK is in- 
ternal to the program. It is accessed by typing 'HELP* at 
any time. This gives the user a full-blown explanation of 
the commands and modes available for each module. 

Ic is necessary to run NASCAP before using TERMTALK. 

A simple object definition will initialize enough files to 
use SUBSET. For instructions on NASCAP use, see the NASCAP 
User's Manual - 1978, NASA CR-159417, August 1978. 

To use TERMTALK on a UNIVAC computer, type 
'@XQT TERMTALK.ABS' , The program will ask for the NASCAP 
file prefix and assign all necessary files. 

On CDG machines, the user must attach NASCAP files 
16, 17, and 21 before running. At the end of each run, file 
number 3 will contain a line printer image of the output. 

SAMPLE RUN 

The following shows a simple TERMTALK example. Note 
that near the end, the mode 'NOTERM' is set. This suppresses 
terminal printing of output data. The user will rely on line 
printer output, without having to wait for terminal printout 
of the HISTORY graph. 


:PXQT TERMTFlLK. ftBS 
WELCOME TO TERMTftLK 


PREFIX PLEftSE ? >MM 
THRNK YOU 


~ notice — THESE FILES WERE i5EhERRTEB UNDER ^OLD*' NHSCfiP 
field RND flux INFDRMftTION WILL NOT BE RVfilLRBLE 

CHOOSE RNY MODULE 

HELP IS RLWRYS RVRILRBLE - TYPE ''HELP^ 

>SIN6LE 


SINGLE CDMMRND OR MODE SET ? 
>15 

CELL NUMBER 15 

‘ENTERED RT -6.5 -3.5 « 


CELL NUMBER 15 

CENTERED RT -6.5 -3.5 

MRTERIRL IS SDLRR 
POTENT I RL * -6.33’?-»‘03 VOLTS 


SINGLE CDMMRND DR MODE SET ? 
>RLSO STRESS 
MODE RESET 

SINGLE CDMMRND DR MODE SET ? 
>RLSO DELTA 
MODE RESET 

SINGLE CDMMRND DR MODE SET ? 
>155 

CELL NUMBER 155 

CENTERED RT -5.5 -S.5 8.0 

MRTERIRL IS GDLDPD 
PDTENTIRL « -6.403+03 VOLTS 
DELTA V s 0.000 VOLTS 
INTERNAL FIELD STRESS * 0.000 

SINGLE CDMMRND OR MODE SET ? 

> LATEST 

LATEST COMMAND DR MODE SET ? 
>LIST 1 45 


VOLTS^METER 


PDTL IN VOLTS 

FDR NRSCAP 

CYCLE 

11... 

TIME = 4.23+03 

SEC 


1 -6.44+03 

2 -6.33+03 

3 

-6.33+03 

4 -6.33+03 

5 

-6 . 33+0 

6 -6.33+03 

7 —6.33+03 

3 

-6.41+03 

3 -6.41+03 

10 

-6.41+0 

11 -6.33+03 

12 -6.33+03 

13 

-6.33+03 

14 -6.33+03 

15 

-6 . 33+0 

16-6.41+03 

17 -6.33+03 

13 

-6.44+03 

13 -6.33+03 

20 

-6 . 33+0 

21 -6.33+03 

22 -6.33+03 

23 

-6 . 33+03 

24 -6.33+03 

25 

-6.41+0 

26 -6.41+03 

27 -6.41+03 

23 

-6 . 33+03 

23 —6.33+03 

30 

-6 . 33+0 

31 -6.33+03 

32 -6.33+03 

33 

-6. 41+03 

34 -6.33+03 

35 

-6.44+0 

36 -6.33+03 

37 -6.33+03 

33 

-6.33+03 

33 -6.33+03 

40 

-6.33+0 

41 -6.33+03 

42 -6.41+03 

43 

-6.41+03 

44 -6.41+03 

45 

—6 . 33+0 


LATEST COMMAND DR MODE SET ? 
•HISTORY 
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HISTORY COMMFtHD DR MDDE SET ? 

>-:3 -4 -5 

POTL IN VOLTS 

TIME : -3 -4 

. — : — 


> 


1 . 0-06: 
4. 0-01 : 
8.3-01: 
1 . 8 + 00 : 
3. 0+01 : 
1.3+02: 
2.3+02: 
1.2+03: 
2.2+03: 
3.2+03: 
4.2+03: 
. PRESS 


:: 1 

-4.66+02 
— 1 . 77+03 
-3. 00+03 
-4.51+03 
-6.60+03 
-7.83+03 
-8.48+03 
-9.23+03 
-9.30+03 
-9.30+03 
-9.30+03 
<CR> 


4*0 

-5. 00+02 
-1.87+03 
-3. 17+03 
-4.74+03 
-6.30+03 
—6 . 26+03 
-6. 17+03 
-5.39+03 
-6.50+03 
-6 . 00+03 
-6. 12+03 


-4.67+02 
-1.75+03 
-2.96+03 
-4.42+03 
-5.63+03 
-4.96+03 
-4.43+03 
-4. 03+03 
-3.71+03 
-3.46+03 
-3. 25+03 


TO CONTINUE 



«5 


4 *^ 


«7 


POTL VERSUS ldg<:time> 

0.00 + 

: 3 

-2. 00+03+ 

-4. 00+03+ 


-6.00+03+ 

■ 

-8. 00+03+ 

-7.00+00 -5.00+00 -3.00+00” -1.00+00 


3 



1 . 00+00 3. 00+00 5. 00+00 


HISTORY CDMMftND DR MODE SET ? 
>CYCLES 6 TO 11 
MODE RESET 


HISTORY COMMPfND DR MODE SET ? 
>GRRPH 
MDDE RESET 

HISTORY CDMMftND DR MDDE SET ? 
>-3 -4 -5 
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PQTL VERSUS LDG<TIME> 
-3. 00-»-OS+ 


3 


-4. 00+0S+ 

-!5. 00+033 

“ 6 . 00403 + 
£ 

-7. 00+034 
-3. 00+031 

-9. 00+03+ 

•- 

9 - 

-1 . 00+04+- 
£. 10+00 


3 


£ 




£ £ 


50+00 


2.90Hh00 


4 »- 

3.30+00 


1 1 
3.70+00 


HISTORY CDflMflHD DR MODE SET ? 

>SUBSET KfiPTON 

DEFIMITION OF NEW SUBSET NRMEl) KftPTOH 
1£67 remaining in GROUP 


? 


SUBSET INSTRUCTION PUEflSE 
>MPITL KfiPTON 

6 REMfilNING IN GROUP 

SUBSET INSTRUCTION PLEfiSE ? 

> WHICH? 

MEMBERS OF GROUP KfiPTON 

90 91 107 103 339 1130 0 0 0 

6 REMfilNING IN GROUP 

SUBSET INSTRUCTION PLEfiSE ? 

>DDNE 

GROUP KfiPTON WITH 6 MEMBERS IS NOW DEFINED 
RETURNING TO MODULE -'HISTDR'’ 

HISTORY CDMMfiND DR MODE SET ? 

>SPECIfiL 

SPECIfiL CDMMfiND DR MODE SET ? 

>NDTERM 
MODE RESET 

SPECIfiL CDMMfiND DR MODE SET ? 

>H I STORY 

HISTORY CDMMfiND DR MODE SET 7 
> GROUP KfiPTON 

HISTORY CDMMfiND DR MODE SET ? 

>EXIT 

WOULD YOU LIKE fi LINE PRINTER COPY OF THIS 
TERMTfiLK SESSION ?>YES 
-EJaT TERMTflLK- 

FIJRPUR £7R3fi E35 73Rl-£ 03/££/79 1G:£4:5£ 

3 BLOCKS COPIED. 

REfiDY 
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INTERNAL DOCUMENTATION 


This is a computer printout of the actual on-line 
documentation available in TERMTALK. First is the complete 
^ output of the OUTLINE aid. Second is a listing of all 

available HELP sequences. 

i 

I 

i 

I 


i 

I 

I 

i 



I 

e: 
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CHOOSE HNY MODULE 

HELP IS flLWPYS flVftlLFIBLE - TYPE ''HELP" 
> OUTLINE 

PICK ft NUMBER FDR THE DESIRED OUTLINE: 


1> 

MODULE 

ftND 

ftID NftMES 


HISTORY 



S> 

LftTEST 




SINGLE 



S’) 

SPECIftL 



G> 

SUBSET 



?> 

ftLL OF 

THE 

ftBOVE 

S) 

NONE OF 

THE 

ftBDVE 

’?> 

CURRENT 

MODULE ONLY 


ftNY ft ID MfiY BE CftLLED FROM ftNY MODULE 

ft IDS 

»»»»»»»♦ 

ftGft I N 
HELP 

LDCftTION «? 

OUTLINE 
SUBSET 
EXIT 

SUBSET CtSRDUP NftMEI 


MODULES 

■»»»»» ♦ »»» 

HISTORY 

LftTEST 

SINGLE 

SPECIftL 


MODES —— 
FLUX - FIELD 
- DELTft 
CYCLE « TO « 
TIME « TO « 
TftBLE - GRftPH 
LINEftR - LOG 


HISTORY MODULE 

COMMftNDS 

POTL « 

*m> I F F f • • V 

GROUP [GROUP NftME] 


- BOTH 
r^UMCYCLE 


LftTEST MODULE 


MODES 

FLUX - FIELD - POTL 

- DELTft - STRESS 
SEQUENT I ftL - MftGNITUDE 
- ftBSMftG 


COMMftNDS 

GROUP [GROUP NftME] 

LIST TO s: 

ftLL 


SINGLE MODULE 


— MODES 

---- 

COMMRNDS 

EVERYTHING - 

NOTHING 

(> 

«:RLSD DR NDT> 

FOLLOWED BYi 

f ONLY COMMRND FDR 

NUMBER 

STRESS 

SINGLE IS TO ENTER R SINGLE 

CENTER 

NDRMRL 

CELL NUMBER > 

MRTL 

SHRPE 


POTL 

CODE 


FLUX 

CNUMB 


FIELD 

CPDTL 


DELTR 




gPECIftL MPPULE 

— — MODES COMMANDS 

TERMPT - NOTERM CYCSET « 


SUBSET SPEC I F I CRT I DNS 


HELP 

INSULRTDP 

BRRE 

DIRECTORY 

WHICH? 

DONE 

OMIT 


XLIM « TO P 
YLIM « TO « 
ZLIM « TO « 
NOPMRL « « « 
BDIREC « a 
NUMBERS TO 
INCLUDE 
CNUMB « 
EXCLUDE 


MRTL CMRTERIRL NRME3 
SHRPE CSCURRE» RECTRN» 

OP CGROUP NRME3 
RND C6RDUP NRME3 
CDMPL CGROUP HRME3 
NRME CRNY WORD3 
CHOOSE RNY MODULE 

HELP IS RLWRYS RVRILRBLE - TYPE '‘HELP 
>HELP 

HELP IS RT HRND - PICK R NUMBER? 

<l> BRSIC USE 

CURRENT MODULE 
NUMBER I N»3 CONVENT IONS 
FLUX? FIELD? POTL? DELTR? 

SUBSET RND GROUPS 
RIDS 

as? RETURN 


RIGHT? ECUIL? OR BDOM3 


<4) 

<5> 

<6> 


STRESS 


TO 


C7? 

m> 

< 9 ? 

< 1 0 ) 
ai> 
as? 

MR I N 


ERROR MESSRGES 
COORDINRTE SYSTEM 
LINE PRINTER FILE 
DEFRULT MO’^ES 
OLD NRSGRP 
CDMPLRINTS 


1 

— ERSIC USE — 

TERMTRLK CONSISTS OF fl SET OF MODULES? ERCH GIVING 
THE USER RCCESS TO R CERTRIN TYPE OF NRSCRP INFDRMRTIDN. 
•‘‘HISTORY-' FDR EXRMPLE? GIVES THE CHRRGING HISTORY OF 
INDIVlnURL SUPFRCE CELLS. ^LRTEST*' GIVES INFDRMRTIDN 
ONLY FROM THE LRTEST CHRRGING CYCLE. •'SINGLE'’ GIVES 


GEDMETRICRL 
CELL. 


RND OTHER INFDRMRTIDN RBDUT RN INDIVIDURL 




WHICHEVEP MOrilJLE VDU fiPE IN? YOU HftVE ftVFl I LftELE ft 
SET OF COMMftHnS ftHD ft SET OF MODES. CDMMftHDS IMITIftTE 
OUTPUT. MODES ftLTER THE FORM OF SUCCEED IMG OUTPUT. 

FIT ftMV TIME IHSTEftD OF TYPING ft COMMftMD 
OR MODE? YOU MftY CHftNGE MODULES OR CftLL FOR ftlD. 
ftLTOGETHER YOU HftVE FOUR CHOICES - TYPE ft COMMftMD r 
ft MODE? THE NftME OF ft DIFFERENT MODULES » OR THE 
NftME OF ftN ft ID. 

ftS ftN EXAMPLE,. YOU ARE IN MODULE SINGLE. YOU 
TYPE THE NUMBER 317. THIS IS ft COMMAND. THE 
TERMINAL PRINTS OUT THE POTENTIAL AND THE SURFACE 
NORMAL OF CELL «t317. YOU WANT MORE. 

YOU TYPE ■'EVERYTHING'’ <ft MODE) AND ON THE NEXT 
LINE TYPE '317'' AGAIN. MORE INFORMATION COMES OUT. 

TO CHftNGE MODULES? YOU TYPE •'HISTORY-'. NOW YOU DON'T 
KNOW WHAT TO DO? SO YOU TYPE ''HELP-'. AND YOU GET THIS 
MESSAGE. 

FOR ANY COMMAND? MODE? MODULE? OR AID? TERMTALK 
RECOGNIZES AN ENTIRE WORD BASED ON THE FIRST THREE 
LETTERS. 

PICK ANOTHER NUMBER -OR- TYPE •'MENU'’ TO SEE YOUR CHOICE 
-OR- <CR> TO RETURN TO MftIN 

-- MftIN — 

IF YOU ARE IN MODULE MftIN? YOU HftVE 
JUST ENTERED TERMTALK. I AM WAITING FOR, YOU TO CHOOSE ft 
MODULE. ALL OF THE AIDS APE ftVftlLABLE ftS WELL. HERE 
IS ft list: 

CHOOSE ANY MODULE NftME ... 

ANY AID MftY BE CALLED FROM ANY MODULE 


MODULES AIDS 

HISTORY AGAIN 

LATEST HELP 

SINGLE LOCATION 

SPECIAL OUTLINE 

SUBSET 
EXIT 

SUBSET [GROUP NAME] 

PICK ANOTHER NUMBER -OR- TYPE •'MENU-' TO SEE YOUR. CHOICE 
-OR- <CR> TO RETURN TO MAIN 


3 

— NUMBER I NG CONVENT I ONS ~ 

TERMTALK HAS PARTICULAR NUMBERING CONVENTIONS FDR 
SURFACE CELLS AND CONDUCTORS. A NASCAP OBJECT IS COMPOSED 
OF N FLAT SURFACE CELLS? M BOOM CELLS? AND L CONDUCTORS. 

N OR M MAY E0IJAL ZERO. 

NEGATIVE NUMBERS < -L THROUGH THROUGH -1 > REFER TO 
CONDUCTORS. THE POSITIVE NUMBERS 1 THROUGH N REFER TO 
FLAT SURFACE CELLS. THESE ARE THE SMALL RECTANGLES AND 
TRIANGLES COVERING THE SPACECRAFT. BOOM CELLS ARE SKINNY 
CYLINDERS. THEY ARE NUMBERED N+l THROUGH N+M. 

... PICK ANOTHER NUMBER -OR- TYPE •'MENU-' TO SEE YOUR CHOICE: 
-OR- '^CR?' TO RETURN TO MAIN 
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4 

— FLUX FIELD PDTL DELTR STRESS — 

THESE RRE THE FIVE DYHRHIC PIECES DF INFDRMRTIDN 
PERTRINIHI3 TO RNY SUPFRCE CELL. THEY CHRNGE ERCH TIME 
CYCLE. OHLY FLUX RHD POTL PERTRIM TO CONDUCTORS. IF 
YOU RRE RHRLVZIMG -'OLD'' NRSCRP RUNS? FLUX RHD FIELD 
WILL HOT BE RVRILRBLE. 


FLUX - TOTRL ELECTRIC FLUX TO R CELL OR CONDUCTOR - 
DUE TO INCIDENT PRRTICLES^ SECONDRRIES-. 

ELECTRON BRCKSCRTTERj RND PHOTOEMISSION. 

UNITS RRE RMPS PER SQURRE METER. 

FIELD - THE EXTERNRL ELECTRIC FIELD IN THE VOLUME 
IMMEDIRTELY RBOVE R SUPFRCE CELL. UNITS OF 
VOLTS PER METER. 

POTL - ELECTRIC POTENT I RL ON THE EXTERNRL 

SURFRCE DF R CELL DR R CONDUCTOR. UNITS 
RRE VOLTS. 

DELTR - DIFFERENCE BETWEEN SURFRCE POTENT IRL DF R 
DIELECTRIC CELL RND POTENT IRL OF THE 
UNDERLYIN6 CONDUCTDR. IN VOLTS. 

STRESS - INTERNRL ELECTRIC FIELD STRESS FOR DIELECTRIC- 
CELLS. .PROPERTIONRL TO DELTR. INVERSELY 
PROPDRTIONRL TO MRTERIRL THICKNESS. UNITS 
DF VOLTS PER METER. 

... PICK RNDTHEP NUMBER -DR- TYPE 'MENU'' TO SEE YOUR CHOICES 
-DR- <CR> TO RETURN TO MRIN 






— SUBSET RND SRDUPS — 

YOU CRN USE SUBSET TO DEFINE R fiRDUP DF CELLS THRT 
RRE DF SPEC IRL INTEREST TO YOU. FDR INSTRNCE-. YOU MIGHT 
WRNT RLL KRPTDN CDRTED SURFRCE CELLS. DR RLL THE CELLS 
OVER CDNDUCTDR DR THE INTERSECTION DF BOTH SETS. 

YOU ENTER SUBSET BY TYPING SUBSET ■- j DR --SUBSET^ 


FOLLOWED BY R GROUP NRME. GROUP NRMES MUST BE TYPED 
OUT CORRECTLY r NOT JUST THE FIRST THREE LETTERS. YOU 
MRY DEFINE UP TO 34 DIFFERENT GROUPS.'. IN ADD IT ION TO 
THE •£ DEFAULT GROUPS -RLL'' RND -''NULL-',. WHICH 
REPRESENT THE COMPLETE SET RND THE EMPTY SET. 

THE SUBSET SPECS CRN EE CLASSIFIED INTO NARROW 
DOWN TYPES? SET OPERATIONS? INFORMATION REQUESTS? RND 
OTHERS. 


NARROW DOWN ***** 

R NEW SUBSET CONSISTS DF RLL CELLS RND CONDUCTORS. 

THE NARROW DOWN SPECIFICATIONS ELIMINATE ANY CELLS THRT 
DON-’T FIT. RLL DF THE NARROW DOWN SPECIFICATIONS? 
EXCEPT •-NUMBERS'’? ELIMINATE RLL CONDUCTORS? -1 
THROUGH -7. 

XLIM « TO « - EXCLUDE CELLS WITH CENTERS OUTSIDE DF 

GIVEN LIMITS IN THE X DIRECTION 
YLIM it TO a - LIKE XLIM 
ZLIM « TO - LIKE XLIM 
INSULATOR - EXCLUDE CONDUCTING CELLS 
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BRRE - EXCLUDE INSULf=lTIM*5 CELLS 

NDRMftL « « « - EXCLUDE RHY CELL WITH R SURFRCB HDRMRL 

□THER THRH THIS. ELIMIMRTE RLL BOOM CELLS. 

BDIREC « « - EXCLUDE BOOMS MOT POINT I MB IN THIS 

DIRECTION. ELIMINRTE RLL FLRT SURFRCE CELLS. 
CNUMB « - EXCLUDE CELLS WITH UNDERLYING CONDUCTORS 

OTHER THRN THIS ONE 

MRTL CMRTERIRL3 - EXCLUDE CELLS NOT OF TH-IS MRTERIRL 
i FIRST 6 LETTERS SI6NIFICRNT > 

SHRPE CSHRPE3 - EXCLUDE RNY CELL NOT IN T>ff ;• SHRPE 
< 3 LETTERS RRE ENOUGH > 

NUMBERS TO « - EXCLUDfa IF CELL NUMBERS RRE NOT 

WITHIN THESE LIMITS. IF THE FIRST NUMBER IS LESS 
THRN -?!* NO CONDUCTORS WILL BE ELIMINRTED. 

SET OPERRTIDNS >»»»♦ 

THESE RLLOW YOU TO COMBINE THE CURRENT GROUP WITH- 
RNY PREVIOUSLY DEFINED GROUP. TO STRRT OVER- ON CURRENT 
GROUP IT '•OP RLL^ OR- '’COMPL NULL-' BRINGS BRCK RLL 
MEMBERS !T WHILE '‘RND NULL'’ DR ''CDMPL RLL'' 

EUMINRTES RLL. THE GROUP NRMED < EXCEPT FDR ''RLL-' 

RND 'NULL'’ ) MUST HRVE BEEN DEFINED BY THE USER. 

DR CNRME3 - PRESENT GROUP BECOMES UNION OF PRESENT 
GROUP RND NRMED GROUP. 

RND CNRME3 - PRESENT (GROUP BECOMES INTERSECTION OF 
PRESENT (GROUP RND NRMED GROUP. 

CDMPL CNRME3 - PRESENT (GROUP BECOMES CDMPLEMEHT OF 
NRMED (GROUP. 

INFDRMRTIDN 

WHICH? - PRINTOUT OF CURRENT GROUP MEMBERS 
DIRECTORY - GIVES LIST OF PREVIOUSLY DEFINED GROUPS 
HELP - THIS MESSRGE 

OTHER SPECS »»■»»» 

NRME CNRME3 - CHRNGE THE NRME OF CURRENT GROUP TO THIS 
INCLUDE [NUMBERS] - UP TO 14 SPECIFIED CELLS RND 
CONDUCTORS BECOME GROUP MEMBERS . 

EXCLUDE [NUMBERS] - UP TO 14 MEMBERS REMOVED FROM' (GROUP 
DONE - CRT(=(LDGUE THIS (GROUP FDR LRTER USE. RETURN- TO 
MODULE. 

OMIT - DON'T CRTRLDGUE THIS (GROUP. RETURN TO MODULE. 


SUBSET 

SPECIFICRTIDNS 

HELP 

XLIM TO 

INSULRTDR 

YLIM o TO « 

BRRE 

ZLIM TO tJ- 

DIRECTORY 

NORMRL 

WHICH’? 

BDIREC 

DONE 

NUMBERS TO « 

OMIT 

INCLUDE .... 


CNUMB » 


EXCLUDE .. . 






MFITL CMftTERIftL NRWEJ 

SHRPE CSQUFtPEi^ RECTRN? RIEHT? EQUILr DR- BODftD 

□R CGRDUP NRME3 

RHD CGRDUP MRME3 

CDMPL C.GROUP HRHE3 

NRME CRNY W0RD3 

... PICK RNDtHER RUIIBER -DR- TYPE •'MENU'' TQ SEE YDUR CHDICES 
-DR- <CR> TD RETURN TD MRIN 

6 

— RIDS — 

RIDS RRE RVRILRBLE FROM RNY MODULE. YDU CRNMDT 
GENERRLLY CRLL DNE RID WHILE YDU RRE USING RNDTHER. 

YOU CRNNOT CRLL ''EXIT^ WHILE USING 'HELP-'r FDR 
EXRMPLE. THE RIDS RRE : 

RGRIN - CRUSES TERMTRLK TD REPERT THE- MDST RECENT PROMPT 
HELP - PRINTS VRRIDUS INFDRMRTIVE MESSRGES 
LDCRTIDN - 'LDC-' DR 'LDCRTIDN^ FOLLOWED BY R 

CELL NUMBER WILL GIVE THE CDDRDINRTES DF THE 
CENTER DF THRT SURFRCE CELL. 

OUTLINE - PRINTS R MENU:. R LIST DF RVRILRBLE 
MODULES-. COMMRNDSj. MODES ^ RND RIDS. 

SUBSET - RLLDWS USER TD DEFINE R GROUP DF CELLS RND 
CONDUCTORS THRT RRE OF PRRTICULRR- INTEREST. 

THE ENTIRE SET DF CELLS IS NRRRDWED DOWN TD 
R DESIRED SUBSET. 

SUBSET CGRDUP NRME3 - IF THE GROUP NRME HRS BEEN 
PREVIOUSLY DEFINED^ PERMITS RLTERRTIDN DF 
THE MEMBERS. IF NDTr R NEW GROUP WITH THRT 
NRME IS DEFINED. 

EXIT - ENDS TERMTRLK USE. 

... PICK RNDTHER NUMBER -OR— TYPE '’MENU'' TD SEE YDUR CHOICES 
-DR- <CR> TD RETURN TD MRIN 

-- ERROR MESSRGES — 

LET'-S SRY YDU DID SOMETHING I DIDN'T EXPECT. OR 
LET'S SRY I JUST GOT CONFUSED. I WILL PRINT OUT RN 
ERROR MESSRGE. THE FIRST WORD DF RN ERROR MESSRGE HRS 
RSTERISKS RROUND IT.* RND IDENTIFIES THE SUBROUTINE WHERE 
THE PROBLEM DCCURED? FDR EXRMPLE SINGLE . 

WHRT DD YDU DO THEN? TRY IT RGRINi. DR TRY SOMETHING 
ELSE. YDU CRN'T DRMRGE THE PROGRRM. DR CRLL FDR THE 
'OUTLINE' DR ^HELP". 

.... PICK RNDTHER NUMBER -OR- TYPE 'MENU' TO SEE YOUR CHOICES 
-OR- <CR> TD RETURN TD MRIN 


8 

— COORDIMRTE SYSTEH — 

THE CODRDIMftTE SYSTEH USED BY TERHTRtK HRS THE OR 1 81 H 
IM THE CENTER OF THE INNERMOST GRID. LOCRTIONS IN RLL 
GRIDS HRE IN TERMS OF INNER GRID UNITS. 

IN OTHER WORDS r PERMITTED VRLUES FOR M DR Y IN THE 
INNER GRID RRE FROM -8 TO ■►S. IN GRID «a VRLUES RRE 
FROM -10 TO -18 IN STEPS OF TWOr DR FROM -*“10 TO -t-l* BY 
TWOS. Z GOES FROM -16 TO -HlG IN THE INNER »5RIDr RND FROM 
-3£ TO -HSa IN GRID «£. 

... PICK RNOTHER NUMBER -DR- TYPE '‘MENU'' TO SEE YOUR CHOICES 
-OR- <CR> TO RETURN TO MR-IN 

9 

— LINE PRINTER FILE — 

TERMTRLK CRERTES R LINE PRINTER FILE. RT THE END 
OF R RUN? YOU CRN SEND THIS FILE TO THE PRINTER TO MRKE 
R HRRD COPY. IF YOU RRE USING R SLOW TERMINRLr YOU CRN 
RSK FDR MODE ^NDTERM" IN MODULE '‘SPECIRL''. IN- TH-IS 
MODE? RLL THE OUTPUT GOES TO THE LINE PRINTER FILEr RND 
NONE DF IT GOES TO YOUR TERMINRL. YDU CRN STILL GET R 
COMPLETE HRRD COPY RT THE END OF THE RUN. 

... PICK RNOTHER NUMBER -OR- TYPE ''MENU'- TO SEE YOUR CHOICES 
-OR- <CR> TO RETURN TO MRIN 

10 

— DEFRIJLTS — 

RT THE BEGINNING DF TERMTRLK RLL MODES RRE SET TO R 
DEFRULT VRLUE. THE DEFRULT MODES FDR ER*:H' MODULE RRE 

HISTORY: PDTL? RLL CYCLES? BDTH? LOG 
LRTEST: PDTL? SEQUENT IRL 
single: RLSD NUMBER? CENTER? MRTL? PDTL 
NOT < EVERYTHING ELSE > 

SPECIRL: TERMPT 

... PICK RNOTHER NUMBER -OR- TYPE 'MENU'' TO SEE YOUR CHOICES 
-DR- <CR> TO RETURN TO MRIN 

11 

— OLD NRSCRP — 

TERMTRLK CRN BE USED ON NRSCRP FILES THRT WERE' 

GENERRTED BEFORE TERMTRLK WRS WRITTEN. BUT IN THE OLD 
DRYS? NRSCRP DIDN'T SRVE FLUX OR FIELD INFORMRTIDN ON 
THE RESTRRT FILES. SO THOSE NUMBERS RREN''T RVRILRBLE 
FROM OLD NRSCRP RUNS. 

... PICK RNOTHER NUMBER -OR- TYPE 'MENU'' TO SEE YOUR CHOICES 
-OR- <CR> TO RETURN TO MRIN 

1£ 

— COMPLRINTS — 

TERMTRLK WRS WRITTEN TO BE R CDNVENIEMT?' USER-ORIENTED 
DRTR REDUCTION TOOL. 

IF YDU DDN^T LIKE IT? DR WOULD LIKE TO SEE SOME 
CHRNGES? YDU CRN WRITE TO SYSTEMS? SCIENCE? RND SDFTWRRE? 

PD BOX 16£0? LR JDLLR? CRLIFORNIR 9£038. MRRK IT 
'RTTENTIDN: JRCK CRSSIDY' 

... PICK RNOTHER NUMBER -OR- TYPE •'MENU' TO SEE YOUR CHOICES 
-OR- <;CR> TO RETURN TO MRIN 
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CHOOSE HHY MDCULE 

HELP IS RLWRYS RVRILRfiLE ~ TYPE '’HELP^ 
>HI STORY 


HISTORY CDHMRHH DR MODE SET ? 

>HELP 

HELP IS RT HRHD - PICK R NUHBERr 
<1> BRSIC USE 
<£> CURREHT MODULE 
< 3 > NUMBERIH6 COHVEHTIDMS 
<4> FLUXfFIELDf PDTLf DELTf^? STRESS 
<5> SUBSET RHD GROUPS 
RIDS 


<7> ERROR MESSRGES 
<8> CODRDINRTE SYSTEM 
<9> LINE PRINTER FILE 
<10> DEFRULT MODES 
<11> OLD NRSCRR 
<1£> CDMPLRINTS 


a3> RETURN TO HISTDR 


£ 

THE CDMMRNDS FDR MODULE ''HISTORY'’ RRE- 
CELL NUMBERS. YOU CRN TYPE ONE CELL NUMBERy DR R STRING 
OF UP TO FIFTEEN NUMBERS ON THE SRME LINE. CONDUCTOR 
NUMBERS <NEGRTIVE> RRE RLSO VRLID. IF YOU HRVE R SET OF 
NUMBERS YOU PLRN TO USE MORE THRN DNCEr YOU CRN DEFINE 
THEM RS R GROUP USING ''SUBSET". 

YOUR OUTPUT WILL BE R TRBLE RND/OR R GRRPH 
OF WHICHESCR ONE OF THE 5 QURNTITIES IS THE CURRENT 
MODE. NOTICE THRT THE "HISTORY" MODES RRE DISTINCT 
FROM THE "LRTEST" MODES. IF YOU SPEICIFY R RRNGE OF 
CYCLES OR R RRNGE OF TIMES» ONLY THOSE ONES -WILL BE 
INCLUDED IN THE TRBLE RND GRRPH. THE X RXIS OF THE 
GRRPH WILL BE THE TIME? THE LOG BRSE 10 OF THE 
TIME? DR THE CYCLE NUMBER - DEPENDING ON WHETHER- 
THE MODE IS "LINERR"? 'LOG''^ DR "NUMCYC". 


HISTORY MODULE 


MODES COMMRNDS — 

FLUX - FIELD - POTL « 

- DELTR - STRESS sjj ... 

CYCLE TO « GROUP [GROUP NRME3 

TIME TO « 

TRBLE - GRRPH - BOTH 
LINERR - LOG - NUMCYCLE 

. .. PICK RNOTHER NUMBER -OR- TYPE "MENU" TO SEE YOUR CHOICES 
-OR- <CR> TO RETURN TO HISTOR 


HISTORY COMMRND OR MODE SET ? 
LRTEST 


LRTEST COMMRND OR MODE SET i 
>HELP 

HELP IS RT HRND - PICK R 
Cl) BRSIC USE 
<£> CURRENT MODULE 
<3) NUMBERING CONVENTIONS 
<4> FLUXir FIELD!. POTL!. DELTRj 
C5> SUBSET RND GROUPS 
<6'> R IDS 

(13) 


UMBERy 

<7) ERROR MESSRGES 
<8> CODRDINRTE SYSTEM 
•::9> LINE PRINTER FILE 
STRESS aO) DEFRIJLT MODES 
ai> OLD NRSCRP 
a£> CDMPLRINTS 
RETURN TO LRTEST 
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LRTEST PRINTS INFOPMRTIDH RBOUT R GROUP OF 
CELLS RND CONDUCTORS. THE THREE COMMRNDS RRE THREE WRYS 
TO SPECIFY R GROUP. '’RLL'’ WILL INCLUDE THE FULL SET 
□F CELLS RND CONDUCTORS. •’ GROUP ^ FOLLOWED BY R GROUP 
NRME PRINTS DNY THOSE IN THE USER-DEFINED GROUP. •'LIST'' 
FOLLOWED BY TWO NUMBERS WILL LIST THE INTERVENING CELLS. 

THE INFORMRTION PRINTED CRN BE FLUXr FIELDf POTLf 
DELTRi* DR STRESS. 

IF THE MODE IS 'SEQUENT I RL^ THE INFDRMGTIDN 
WILL BE ORDERED BY CELL NUMBER. IF •' MRGN I TUBE ' f IT 
WILL BE ORDERED BY THE VRLUES THRT RRE BEING PRINTED. 

FDR MODE 'RBSMRG^f ORDERING IS BY RBSDLUTE VRLUES. 


LGTEST MODULE 


MODES 

FLUX - FIELD - PDTL 

- DELTR - STRESS- 
SEQUENT I RL - MRGNITUDE 
- RBSMRG 

. . PICK RNDTHER NUMBER 
-DR- <CR> TO RETURN 


COMMRNDS 


GROUP CGRDUR 
LIST » TO « 
RLL 


-DR- TYPE 
TO LRTEST 


'MENU'' TO 


NRNE3 


SEE YOUR CHOICE 


LRTEST CDMMRND DR MODE SET 1 
> SINGLE 

SINGLE CDMMRND DR MODE SET ? 

>HELP 

HELP IS RT HRND - PICK R NUMBER » 


a> BRSIC USE 

<2> CURRENT MODULE 

<3> NUMBER I NG CONVENT I DNS 

<4> FLUX ..FIELDf PDTLf DELTRr STRESS 


•:5> SUBSET RND GROUPS 
C6> RIDS 

RETURN- TO 


-- SINGLE — 


<7> ERROR MESSFH5ES 
<8> CDDRDINRTE SYSTEM 
<-rK) LINE PRINTER FILE 
ao> DEFRULT MODES 
<11> OLD NFkSCRP 
<1S> CDMPLR-INTS 
SINGLE 


THIS MODULE GIVES VRRIDUS INFORMRTION RBDUT R 
SINGLE FLRT CELL DR BOOM CELL. THE ONLY CDMMRND IS THE 
NUMBER OF R CELL. 

MODES SPECIFY WHICH- CELL PROPERTIES RRE TO BE PRINTED. 
IF YOU WRNT RLL DF THEMf SET 'EVERYTHING-'. TO WIPE THE 
SLRTE CLERNf TYPE 'NOTHING'. IF YOU WRNT TD INCLUDE 
R CELL PROPERTY f TYPE 'RLSD' FOLLOWED ON THE SRME LINE 
BY THE PROPERTY NRME - FDR EXRMPLE 'RLSD POTL'. TD 
EXCLUDE R PROPERTY f TYPE 'NOT' RND THE PROPERTY NRME. 

THE DYNRM I C PROPERTIES f WHICH CHRNGE FROM ONE CYCLE 
TD THE NEXTf RRE 'FLUX' r 'FIELD' r -POTL'r 'DELTFKf 
RND 'STRESS'. 




THE STfiTIC PROPERTIES RRE: 

NUMBER - THE CELL NUMBER 

CENTER - CODRDINRTES OF THE CELL CENTER. FOR BOOM 
CELLS THIS IS MIDWRV BETWEEN THE ENDPOINTS. 

MRTL ~ MRTERIRL NRME 

NORMRL - DIRECTION OF THE NORMRL POINTING OUT OF R FLRT 
SURFRCE CELL. FOR BOOMS » PRINTS THE DIRECTION 
Qp BDDM 

SHRPE - CELL SHRPE. "eITHER 'SQURRE''> 'RIGHT TRIRNGLE"? 

RECTRNGLE '' y '' EQU I LRTERRL TR I RNGLE ? OR ^ BOOM . 
CODE - INTERNRL CELL CODE USED BY NRSCRP. R 36 BIT 
OCTRL NUMBER. 

CNUMB - CONDUCTOR NUMBER OF UNDERLYING CONDUCTOR. 

CPDTL - POTENTIAL OF UNDERLYING CONDUCTOR. 


SINGLE MODULE 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

MODES CDMMRNDS 

EVERYTHING ” NOTHING « 

<RLSO OR NOT!) FOLLOWED BY: < ONLY COMMAND FOR 

NUMBER STRESS SINGLE IS TO ENTER A SINGLE 

CENTER NORMAL CELL NUMBER > 

MATL SHRPE 

POTL CODE 

FLUX CNUMB 

FIELD CPOTL 

DELTA 

... PICK ANOTHER NUMBER -OR- TYPE '^MENU^ TO SEE YOUR CHOICES 
-OR- <CR> TO RETURN TO SINGLE 


SINGLE COMMAND DR MODE SET ? 

*> ;*;PEn I PL 

’ SPECIAL COMMAND OR MODE SET ? 


>HELP 

HELP IS AT HAND - PICK R 

a;> BRsrc use 

< 2 > CURRENT MODULE 
<3> NUMBER I NG CONVENT I DNS 
< 4 > FLUX? FI ELD? POTL» DELTA !• 
<5> SUBSET AND GROUPS 
<6> RIDS 

<13> 


lUMBERf 

<7> ERROR MESSAGES 
<8> COORDINATE SYSTEM 
< 9 > LINE PRINTER FILE 
STRESS <10> DEFAULT MODES 
<n:) OLD NASCAP 
< 12 :> COMPLAINTS 
ETURN TO SPECIA 
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p 


— SPECIRL ~ 

THIS COULD BE CRLLED MODULE MISCELLRMEOUS. 
PRBSEMTLV^ IT HRS ONLY TWO FUHCTIONS - PRINT 
CONTROL RND CYCLE RESETTING. 

RLL PRINTOUT GOES BOTH TO THE TERMINRL RND TO FILE 
3. IF YOU RRE GENGRRTING R LOT OF PRINTOUT RND DDN^T 
WRNT TO WRIT FDR IT RT THE TERMINRL^ USE MODE '‘NDTERM'’. 
THEN YOUR PRINTOUT WILL GO ONLY TO FILE 3r WHICH- CRN BE 
SENT TO THE LINE PRINTER RT THE END OF TERNTRLK'. 

YOU USE CYCLE RESET IF YOU RRE GOING TO MRKE 
FURTHER NRSCRP RESTRRT RUNSr BUT HRVE NO FURTHER USE 
FDR THE DRTR RLPERDY GENERRTED. IN •'SPECIRL" IF 
YOU TYPE '•CYCSET RLL OLD DRTR WILL BE DISCRRDED 
FROM THE •'HISTORY'- FILES. ^LRTEST'* RND '‘SIN»5LE'‘ 

WILL REMRIN UNCHRNGED. 


SPECIRL MODULE 


MODES COMMRNDS 

TERMPT - NOTERM CYCSET o 

... PICK RNOTHER NUMBER -OR- TYPE '’MENU'’ TO SEE YOUR CHOICES 
-OP- 'iCR> TO RETURN TO SPECIR 

SPECIRL COMMRND OR MODE SET ? 

>EXIT 
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5.2 


MATCHG 


Bulk oonductivity currents and a facility to specify 
completely aligned fluxes have been added to MATCHG. Also/ 
MATCHG now handles multiple maxwellian plasmas. These 
changes and two minor input-output changes are described 
here. 

5.2.1 Bulk Conductivity 

For insulating materials, MATCHG now includes the 
bulk conductivity current, J^on' calculation of the 

net flux: 

^oon = - (l) (v - - CE (5.1 

where a is the bulk conductivity, d the material thickness, 

V the surface potential, and the fixed potential of the 
conducting back plate underlying the insulator. MATCHG will 
now request the user to specify for runs with insulating 
materials. The bulk conductivity current will be significant 
whenever (a/d) is of comparable magnitude to ' where 

Jnet incident flux from the plasma or gun. Note 

that equilibrium potentials in MATCHG can now vary as the 
material thicknesses are changed. This new feature of 
MATCHG will allow a crude prediction of insulating surface 
potentials as a function of satellite ground potential. 

5.2.2 Aligned Fluxes 

MATCHG will now allow simulation of either an isotropic 
or a completely aligned incident particle distribution. For 
the isotropic case , the material collects current as would a 
part of a spherical probe. For the completely aligned case, 
the user will be asked to specify 0, the angle of incidence 
with respect to the material surface normal (0® <0 < 90®; 
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default is 6 ■ 0®). The material will then collect current 
as could part of an infinite plane, i.e., for attracted 
particles, 

•^inc * <5-2) 

where is the incident proton or electron flux, T the 

temperature, and N, e, and m the particle density, charge, 
and mass respectively. Unlike the isotropic case, the 
planar collection model has no dependence on the surface 
potential for attracted particles. For repelled particles, 

Jinc is reduced by a factor as in the spherical 

case. 

In the tank charging case, the flux is assumed com- 
pletely aligned, and MATCHG will now always request the user 
to specify the angle of incidence, 9. The incident flux will 
then be = Jq cos0 , where is the flux from the gun at 

normal incidence. 

The formulation used to calculate electron backscatter 
in MATCHG is chosen to be consistent with the incident flux 
direction; an average over a cos0 distribution is used for 
isotropic fluxes, and the yield at a fixed angle 0 is cal- 
culated for aligned fluxes. If the 'ANGLE* secondary formula- 
tion is used, the above statement applies to proton and elec- 
tron secondary yields as well. Use of the 'NORMAL* secondary 
formulation overrides these implicit choices and forces the 
secondary yield at normal incidence to be used irrespective 
of the incident flux direction. Note that the use of 'NORMAL* 
is somewhat inconsistent with an aligned flux at 0 0°; 

MATCHG will print a warning when this choice is made. 

In summary, specification of the MATCHG flux formula- 
tion requires choice of incident flux type, incident flux 
direction, and s^soondary emission formulation, as indicated 
here : 




A. INCIDENT 
FLUX 
TYPE 


TANK 


B. INCIDENT 
FLUX 

DIRECTION 

ISOTROPIC 


C. SECONDARY 
EMISSION 
FORMULATION 

ANGLE 


MAX 

(single Maxwellian) 


ALIGNED 


NORMAL 


DBLMAX 

(double Maxwellian) 


All combinations of the descriptors in columns A^ B/ and C 
are allowed except that the TANK flux type is assumed to be 
aligned. 


5.2.3 Emission Table 

The table of selected emission values produced by 
MATCHG is now consistent with the emission formulations used 
in charging. For a given energy in the table, MATCHG cal- 
culates the yield at a fixed angle or averaged over a cos0 
distribution, as described above. Previously, the yield at 
normal incidence had been given in the table irrespective 
of the formulation used during charging. 

5.2.4 Initial Voltage 

The latest version of MATCHG requests both the initial 
surface potential and the backing plate potential in keV, not 
eV. 


5 • 2 . 5 Multiple Maxwellian Plasmas 

MATCHG has been revised to accept ambient Maxwellian 
plasmas with up to five distinct electron and ion components. 
The new mode is activated by the keyword 'MLTMAX' during flux 
specification; MATCHG will then prompt the user for the 
number of components. Using 'MLTMAX* with one or two com- 
ponents is equivalent to the old 'MAX' or 'DBLMAX', 


respectively. (The old keywords 'MMC' and ’DELMAX' have been 
retained , however , ) 


There are no default values for the 
components of the Maxwellian description, 
the first two components are now n^ » n^^ » 
1 keV, T ^2 * 


third through fifth 
The defaults for 


__“3 
cm f 


T , » 
el 


- 


5.3 WORKSHOP CODES 

In preparation for the NASCAP workshops held at NASA/ 
LsRC, some special new codes were developed. These make use 
of some Univac features to simplify NASCAP running for the 
new user. The three codes are called FILES, OBJCHECK, and 
POTPLT . 

To run NASCAP on a Univac, a user needs three input 
files (OPTIONS, OBJECT, and FLUX), as well as six scratch 
and restart files numbered 10, IS, 16, 17, 21, and 27. 

PILES will generate executive commands to assign, delete, or 
copy a complete set of these files. 

OBJCHECK generates a three-dimensional picture of a 
satellite or test object. It uses the NASCAP shadowing 
routines, without invoking the rest of NASCAP. It is a 
diagnostic tool at object definition time. 

POTPLT is a NASCAP post processor that produces poten- 
tial contour plots. The user specifies for each plot the 
number of grids desired and the location of the plane where 
the contours are drawn. 

5.3.1 FILES 

FILES will create or destroy a set of NASCAP files, 
or copy from one set into another set. A set of files is 
identified by a file prefix, from one to nine characters 
long. For example, if the user chose 'SCATHA' as a prefix. 


the files would be named ’SCATHAIO’, 'SCATHA15S 'SCATHA21', 
and so on. 

FILES is actually a set of three programs ~ 

FILES. ASSIGN (creation ) , FILES. delete (destruction ) , and 
FILES. COPY (copying) . 

The program asks for the prefix name, and it asks 
whether the files are to be public, private, or temporary. 
Public files are default. A private file can be referenced 
only by the person who created it. Temporary files disappear 
when you sign off the computer. 


WELCOME TO ASSIGN ... THIS PROGRAPI WILL ASSIGN 
A SET OF HEU RESTART FILES FOR NASCAR 
I NEED A FILE PREFIX/ 1=9 CHARACTERS LONG/ NO BLANKS 


PREFIX PLEASE ? 
j j c 


THANK YOU 


PUBLIC OR PRIVATE OR TEMPORARY? 
Pub I i 0 

OASG/ PU UJC-fO . 

QASG.PU JJC1S . 

QASG/PU JJC^S////-iOOO . 

@ASG/PU JJC'iF . 

QASGYPLI JJC2i . 

GASG/PU JJC27 . 

QASG/PU JJCOBJ . 

CaASG/PU JJCFLX . 

GASG/PU uJCOPT . 


ALL 9 FILES SUCCESSFULLY ASSIGNED 
Figure 5.2. An example of use of the FILES. ASSIGN program. 
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5.3.2 OB JCHECK 




OB JCHECK allows the user to create an object and look 
at a 3-D picture of it, before worrying about running NASCAP 
itself. OBJCHECK produces material plots, building block 
plots, and hidden line surface cell plots. 

Before defining the object, OBJCHECK asks for NZ, NG, 
XMESH, and the file prefix. NZ is the number of grid units 
in the Z direction. 33 is standard, but for a smaller grid, 
29, 25, 21, or 17 are allowable inputs. NG is the number of 
nested grids necessary to hold the object, including booms. 
XMESH is the NASCAP zone size. 

OBJCHECK then reads the '-OBJ.* file, where repre- 
sents the user supplied prefix. It echoes all input. 

Finally it asks whether object illustration plots are 
desired. If so, it asks if you want material plots. Mate- 
rial plots are fairly time consuming, and are not needed in 
the early stages of object definition. The last input is the 
direction-of-view for the 3-D object illustration plots. 

The directions-of-view are input as (X,Y,Z) vectors, 
not necessarily normalized. Permutations of plus or minus 
1, 2, and 3 usually give good pictures. After you have as 
many views as you want, a carriage return will end the pro- 
gram. As a final task, OBJCHECK runs a few more tests on 
the object, having to do with cell connectivity. 


n 




Qx«it objcheok. 
ASGFIL; FILE 21 
ASGFIL: FILE i7 


ASGFIL: FILE 3 
NZ> 

33 

NG> 

2 

XMESH > 


ISTAT= 100000000000 
ISTAT= 100000000000 
ISTATs 100000000000 


PREFIX PLEASE ? 
j j c 


THANK YOU 


ADDITIONAL OUTPUT UILL BE ON UNIT 3 
Input echoing omitted from figure. 

Figure 5.3. 


Example of OBJCHECK use. Some printout has 
been omitted for brevity. Object definition 
read from JJCOBJ file. 
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CALLING GENMTL 

) INSLST -- 219 INSULATING SURFACE CELLS FOUND 

234 INCLUDING BOOM CELLS 


SOLAR HAS SURFACE RESISTIVITY OF 
WHITEN HAS SURFACE RESISTIVITY OF 
KAPTON HAS SURFACE RESISTIVITY OF 
ASTROO HAS SURFACE RESISTIVITY OF 
TEFLON HAS SURFACE RESISTIVITY OF 
BOOMAT HAS SURFACE RESISTIVITY OF 


1.0+19 OHMS 
1.0+13 OHMS 
1.0+16 OHMS 
1.0+11 OHMS 
1,0+16 OHMS 
1.0+11 OHMS 


FNDSCE -- 493 SURFACE CONDUCTING EDGES FOUND 

BOOMEJ — S26 EDGES FOUND 

1116 ENTRIES IN REVISED VTXL 
NC DETERMINED BY GETNG TO BE 4 

END GENMTL 

cexit: 

Figure 5.3. Concluded. Ghecking cell connectivity. 
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5.3.3 POTPLT 


POTPLT produces the same type of potential contour as 
NASCAP. The difference is that with NASCAP, you have to de- 
cide before each timestep which contour views you want to 
see. With POTPLT, you can choose a view, look at it, then 
go back and produce some more views, without running NASCAP 
again . 

The user can have the plane of view perpendicular to 
the X, Y, or Z axis, located at any node point. Each picture 
can include anywhere from one grid up to the complete computa- 
tional space. 

POTPLT first asks for the file prefix. It then gives 
the options of smoothed contour lines (as opposed to quick, 
rougher lines) and of absolute, 1 to 33, coordinates (as 
opposed to centered, -16 to +16 coordinates) . A carriage 
return will bypass these. 

For each picture the user specifies number of grids, 
view direction, cut value, and number of contours. The 
number of grids can be from I to NG, the total number of com- 
putational grids. The view direction is X, Y, or Z, and is 
perpendicular to the contour plane. So if Y is the view 
direction, contours are drawn in the X-Z plane. The cut 
value determines where the contour plane will cut through 
computational space. The number of contours is set to a 
default of approximately 20. You may choose more or less. 

After all desired views are generated, respond to 
the 'GRIDS TO PLOT <EXIT> ' question with a carriage return 
to end the program. 
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HI? 


PREFIX PLEASE ? 
j J c 


THANK YOU 


TYPE ^SMOOTH' FOR SLOW DRAWING/ SMOOTH CONTOURS 
TYPE 'POSITIVE' FOR POSITIVE COORDINATES Ci-SSI 
<CR> FOR DEFAULT - ROUGH AND CENTERED 

FAST CONTOURS WILL BE DRAWN 
CENTERED COORDINATES C-i6 TO +-163 
GRIDS TO PLOT CEXIT0> 

DIRECTION > 
z 

CUT VALUE C03 > 

3 

HOW MANY CONTOURS? C203 > 

GRIDS TO PLOT CEXITO) 

2 

DIRECTION) 

•< 

CUT VALUE 003 > 

-4 

HOW MANY CONTOURS? 11203 > 

'iO 

GRIDS TO PLOT CEXIT3) 


YOU HAVE created’ 4 FRAMES OF MICROFILM 

CEXIT3 


Figure 5.4. Example of POTPLT use. User generated two con- 
tour planes, at Z = 3 and at X = -4. 


6. PRELIMINARY LOW EARTH ORBIT (LEO) MODEL 

This section describes development of a preliminary 
version of NASCAP designed to simulate charging response of 
high voltage satellites in low earth orbit. During magneto- 
spheric substorms at GEO, typical plasma densities and tem- 

- 3 

peratures are 1 cm ' and 10 keV, respectively, so that the 
plasma Debye length is on the order of several hundred meters. 

For satellites whose characteristic dimensions are a few 
meters, current collection is orbit limited, and spherical 
probe formulations can be used. Furthermore, the charge 
stored on the satellite surface is much larger than the am- 
bient space charge, so that plasma screening can be neglected 
or approximated by linear screening expressions. The original 
NASCAP code made use Of these approximations to simulate 
charging response at GEO. In LEO, the plasma is much denser, 

typically 10 cm , and characterized by a lower temperature, | 

I 

say 1 eV. Debye lengths are on the order of 1 cm, and the | 

orbit limited current collection approximations and linear 
screening expressions are invalid. Therefore, the NASCAP/LEO 
code required the development of sheath limited current col- 
lection algorithms as well as solutions to the Poisson prob- 
lem for the case of highly nonlinear ambient screening. 

A detailed description of the physical models and com- 
putational techniques of the NASCAP/LEO code is contained in 
a paper presented at the AIAA Aerospace Sciences Meeting, | 

January 1980, included in this report as Section 6.4. This 
paper includes a simulation of the experiments performed by 
Konradi and McCoy at the Johnson Space Center. The re- 
sults indicate that ion focusing Observed in the laboratory 
during high voltage collection experiments is probably due 
to voltage gradients on the collecting surfaces. 



Very large ('vl 00-1000 m) spacecraft which have high 
voltage solar arrays are now in the design phase. NASCAP/LEO 
includes a fine resolution mode which allows the current col- 
lection on a periodically continued voltage pattern for a 
large object to be treated with high resolution. This mode 
has been used to study the effects of solar voltage array 
patterns on plasma power losses; our preliminary calculations 
suggest that plasma power loss should not be a primary con- 
sideration in designing the physical arrangement of high 
voltage arrays. These results have been published in IEEE 
Transactions on Nuclear Science, Volume NS-ii7, No. 6, Decem- 
ber 1930, and the paper is Reprinted in Section 6.3. 

Section 6.1 describes calculations comparing NASCAP/LEO 
results to Langmuir probe theory as a simple validation of 
the code. NASCAP/LEO is currently in a very preliminary form, 
and many of the features available in the original NASCAP code 
have yet to be implemented. Section 6.2 assesses the status 
of the code by describing the principle restrictions. 

6.1 LANGMUIR PROBE THEORY COMPARISON 

As a validation of the NASCAP/LEO code, current col- 
lection to a simple object has been calculated for various 
values of the object potential and plasma Debye length. The 

results are compared to the exact numerical calculations for 
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a spherical probe as given by Laframboise. ^ 

The object used in the calculations was a cube with 
1.8' m edges; the zone size was 0.3 m. Plasma temperature was 
1 eV in all cases, and the sheath boundary was defined to be 
at a potential of 0.1 volt. The object voltage was 10 or 20 
volts, and the Debye length was varied from 1 to 100 cm. The 
results are given in Table 6.1. 












TABLE 6.1. LANGMUIR PROBE THEORY COMPARISON 


Ajj (cm) 

V = 

10 Volts 

V * 

20 Volts 

100 

15.7 

(9.6) 

19.7 

(17.1) 

10 

3.1 

(3.8) 

3.8 

(4.9) 

1 

2.7 

(1.9) 

3.1 

(2.0) 

Ratio of 

collected 

current to 

thermal 

current 


culated using NASCAP/LEO for various values of Aq 
(D ebye length) and probe voltage. Numbers in paren- 
theses are for a spherical probe as calculated in 
Reference 7. 

The agreement between NASCAP/LEO calculated values and 
the spherical probe results is good. The effective radius of 
the cubical object employed in the calculations is 1.24 m. 

The ratio of object radius to electron Debye length therefore 
ranges from 1.24 to 124 in these calculations, covering the 
entire range explored by Laframboise. The results are some- 
what sensitive to the potential value used to define the 
sheath boundary. Collected current varies logarithmically 
with the sheath boundary potential; a decade change in the 
boundary potential changes the current by ±50 percent. The 
agreement with the probe theory results is within 50 percent 
in all cases, even though the code resolution ‘is three times 
the Debye length in the = 10 cm case. The cases considered 
for comparison here were chosen to match the parameter ranges 

r 7 1 

reported by Laframboise. ^ The agreement would certainly be 
better if the comparison had been made for a higher object 
potential, whete the physical models in NASCAP/LEO are more 
nearly exact* Given the uncertainty associated with compari- 
son of calculations for a cubical object to a perfect sphere, 
and considering the large zone size, these results provide an 
excellent validation of the nonlinear screening and current 
collection algorithms in NASCAP/LEO. 


6.2 NASCAP/LEO LIMITATIONS 


The present version of the NASCAP/LEO code is intended 
as a preliminary version only. Considerable development will 
be required to bring the code to a state of development similar 
to that of the present NASCAP code itself. The calculations 
described represent the full range of capabilities of the 
NASCAP/LEO code. 

The limitations and restrictions on the preliminary 
NASCAP/LEO model include the following: 

Object Definition . Only objects composed of complete cubes 
are allowed. Booms , flat plates, wedges, tetrahedra, and 
truncated cubes are not provided for. 

Material Properties : Provision for material parameters has 

not yet been included. Secondary emission, photoemission, 
conductivity effects, etc. are therefore not implemented. 

This restriction effectively limits objects to those con- 
sisting of surfaces at fixed potentials. 

Time Dependence ; Only the equilibrium solution is calculated. 
Time dependent charging effects are not simulated. 

Object Potentials ; Floating potentials are not implemented. 

Code Flexibility ; The preliminary version of NASCAP/LEO is 
not yet a user-oriented code in the style of NASCAP. The 
code presently exists in three distinct sections: object 

definition, potential solution, and current collection. 

These sections are executed sequentially, and the user must 
monitor the progress of the results to be sure that proper 
convergence is attained. 


m 







6.3 THE EFFECT OF SOLAR ARRAY VOLTAGE PATTERNS ON 
PLASMA POWER LOSSES 

This is the paper published in IEEE Transactions 
Nuclear Science ^ Volume NS-27« No. 6f December 1980. 




lEEG Tfaiisactions on Nuclear Science) Vol. NS'27, No. 6, December 1980 

THE EFFECT OF SOLAR ARRAV VOLTAGE PATTERNS 
ON PLASMA POWER LOSSES'* 


M. J. MandeXl, Xi Kats> P< G. Sboen, G. W. Schnuelle 
Systems > Science and Software 
Pi 0. Box 1620, La Jolla, CA 92038 


Introduction 


There has been considerable interest re- 
cently in using high voltage solar arrays in 
space, The primary reason is that as power 
requirements increase, cables and PC to DC 
convertors become significant weight factors 
in satellite design, A principal difficulty 
in using high voltages in Low Earth Orbit 
(LEO) , or with a Solar Electric Propulsion 
Stage (SEPS), is the draining of array power 
by currents flowing between exposed surfaces 
through the surrounding plasma, This plasma 
may be of natural origin (as in LEO) , or gen- 
erated by the spacecraft (for SEPS) . Early 
studios by Kennerudi have shown that small 
pinholes in insulation can make almost per- 
fectly insulated solar arrays collect current 
from the surrounding plasma as effectively as 
bare metal plates. Therefore, simply insu- 
lating exposed high voltage areas may not 
prevent substantial power losses. 

Experimental and theoretical studies 
have analysed plasma-array interactions for 
cases in which the potential variation was on 
the scale of the dimensions of the test object, 
For some cases f the magnitude of the power 
loss would prohibit the system from function- 
ing. It has been suggested^ however^ that 
the power loss might be reduced by arranging 
solar cell strings in repeated small-area 
modules to eliminate any large areas at high 
potentials, since the sheath height and area 
increase with the surface potential, smaller 
sheath voltages might reduce total power 
losses. In this paper, the effect of this 
suggestion is examined both analytically and 
with the use of a new, nonlinear plasma sheath 
simulation code. The results indicate that 
small repeated voltage p,ttterns are of mar- 
ginal utility in reducing power loss to the 
surrounding plasma. 

T heory 

The analysis presented is valid for 
plasma current collection by surfaces at high 
voltages. We define high voltages, to be 
voltages such that the associated electron 
potentials, e4>f are much greater than their 
thermal kinetic energy, 

» M f (1) 

There are three length scales of concern; the 
object size, L, the characteristic length for 
potential variation, 1, and the plasma Debye 
length, We require that \ be less than the 
typical object size, Lt 

X < h . ( 2 ) 

This clearly places us in the space-charge- 
limited regime. 

*This3 work supported by NASA/Lewis Research 
center, Cleveland, OR, and Air Force Geo- 
physics Laboratory, Hanscom AFB, MA, 

Contract NAS3-2I762. 


In the limit of very small Debye length 
and very high voltages, a quasi-one-dimen- 
sional sheath forms in front of the collecting 
surface. The current through this sheath is 
the one-sided thermal current, 

" "o 

To find the current structure within the 
sheath requires calculating particle trajec- 
tories and solving Poisson Vs equation in a 
mutually consistent fashion. It is far more 
convenient to approximate the space charge 
density, o, as a function of potential, 
only. The most common such approximation 
(linear, Debye screening), valid for <(> < 0, is 

P(<>) - -e X'' - -2^ . (4) 

ne 

Another useful expression (thin sheath approx- 
imation) is valid for fp >> 0 and locally one- 
dimensional geometry. Requiring the current 
to remain constant at the value given by 
Equation (3) yields 



This expression leads to a nonlinear Poisson 
equation, as discussed below. 

The power loss for an object whose poten- 
tial varies from 0 to is given by 

ma A 

p . (nJAl (vv„^) (6) 

where A is the area of the sheath boundary, 
n is the fraction of charged particles enter- 
ing the sheath that strike the object, and 
^'^max the mean potential at which they 

strike. We would expect v to depend primarily 
on % and X, while A and n depend primarily on 
L and X. If the screening is nonlinear, A 
and V will both depend on The para- 

meter n will usually be near unity. 

Calculations 

The complete self-consistent set of elec- 
tron kinetic equations and Poisson ^3 equation 
are formidable to solve in all but the sim- 
plest cases. Instead we examine here two 
different simplifications which enable us to 
analyze the dependence of the mean collection 
potential on system parameters. 

Harmonic Potential with Linear Screening 

A readily solved case is a large plane 
with potential 

\W 
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V{x,0> ■ y V^JJ a + SinKX) (7) 


wh«r« <L >> 1. If w« rxplAc* Poisaon's aqua- 
tion with a Halmholtz aquation 

(7^ - A“^) V(x,z) - 0 , (8) 

wa find 

V(x#z) ■ y ^ainicx a ^ + a”*'^^j (9) 

2 2 -2 

whara ■ « + X , Tha function v(tfX), 

obtained by numerically calculating particle 
trajactorica, ia ahown in Figure 1. Thia 
function haa a maximum for <X v 1. The total 
variation in v ia from 'v.SO for amall X to 
'v.fS at maximum. Thia would indicate that 
variationa in panel daaign raault in at moat 
a 30 percent variation in power loaa. For 
kX >> 1, particle defleotiona are amall 
becauae the eranaverae fielda are encountered 
at high energy. For i«X << 1/ particle deflec- 
tiona vaniah with the aheath thiekneaa. 


Thraa-Dimenaional Array with Nonlinear 
Screening 


To make a more realiatic eatimata of tha 
plaama loaa. a batter formulation of the 
plaama acreening ia neceaaary aa well aa a 
thrae-dimgnaional treatment of the problem. 

For caaea in which >> 9. a better approx- 

imation for Poiaaon*a equation than the 
Helmholtz equation ( 8 ) ia^ 




7 


1 + (e9/<9)^^^ 


-1 


( 10 ) 


Thia formulation reducea to Equation (4) for 
low potentiala and to Equation (5) for high 
potentiala. The Low Barth Orbit (LEO) cede^ 
ia deaigned to calculate current eoHtiction 
uaing space potentials obtained by solving 
Equation (10). 


The LEO code also haa an expanded resolution 
mode which may be used to calculate v for a 
finally-resolved periodic surface region. 

Zn our previous work^ we calculated the 
sheath structure and power loss for a simu- 
lated Solar array intended to represent t)ie 

experiments of McCoy and Konradi.^ The model 
(Figure 2) consisted of a panel 8.7 m long 
and 1.9 m wide* with a 1.2 m wide high voltage | 

current collecting atrip along its length. I 

Of the several cases previously considered. | 

We will focus on the "linear.” high-voltage \ 

case, in which the strip was linearly biased 
4890 volts from end to end. and on the com- 
parable "uniform” case in which the strip was 
uniformly 2400 volts. The results are shown 
in Table Z for plasma temperature 9 ■ 1 eV 

and Debye length X ■ lo"^ or lo"^ meters. iP„e 
power loss haa been normalized to the plaama 

thermal current («10“^ or lO"^ A/n‘* for an 
oxygen plasma with the respective Debye 
lengths) . 

We wish to compare the power loss from 
the "linear” array with that from a "modular" 
array in which the 0 - 4800 volt potentials 
are contained in 30-cm-square modules (Figure 
3) . Zt is assumed that the gross structure 
of the sheath is the same as that for the 
"uniform" array at the average potential. To 
determine the change in sheath structure near 
the surface, wa used the expanded resolution 
(.0137S m zones) periodic veriion of the LEO 
code. This was done for two different module 
patterns (Figure 3) . The results are given 
in Table ZZ. The power losses are calculated 
by first doubling ;*>e results of cases ZZZ 
and XV of Table I (since increases from 

2400 volts to 4800 volts) . and then multiply- 
ing by '^expanded' resultant losses are 

virtually identical to tha linear cases. Mote 
that, when t.he bracketed factor of Equation 
(10) is taken into account, we are. for these 
parameters, in the region xX > 1 of Figure 1. 

Thus making the pattern finer will reduce 

'^expanded 



Piguro 1. Effective voltage factor, v), for 
sine wave potential [Equation 7) 
as a function of screening length * 



Figure 2 . "Linear” array test object, indi- 
cating bias potentials. The 
labeled portions represent a con- 
ductive strip, while the border, 
sides, and back are insulating 
plastic. For the "uniform” test 
object, the entire conductive strip 
is at 2400 volts, 
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Tftblo I. Powor ho$a 


case 

Potential 

Urn) 

2 

A(m^) 

V 

2 

kw-m Vamp 

I 

0-4800V 

(linear) 

0.01 

1,0 

80.7 

,6935 

269 

n 

0-4800V 

(linear) 

O.OOl 

1,0 

44.4 

.6675 

142 

III 

2400V 

(uniform) 

0.01 

l.O 

92.1 

,975g 

216 

IV 

2400V 

(uniform) 

0.001 

1,0 

48.3 

.9739 

113 





0 Volte 
800 Volte 
2000 Volts 
3200 Volts 
4800 Volts 



Pattern A - "Checkerboard" Pattern B - "Siirip" 

Figure 3, Voltage patterns for module high voltage array power loss calcu- 
lations. Each pattern represents a 30-cm sguare module. 


Table II. Power Loss Parameters - Expanded 
Resolution for .Modular Casas 


Pattern \ (m) Expanded P/J (kw-m^/amp)'*' 


A 

(Checkerboard) 

0.01 

.6210 

268 

A 

(Checkerboard) 

0.001 

.6278 

XA2 

B 

(Strip) 

0 . 01 

,6052 

261 

B 

(Strip) 

0.001 

.6464 

146 


’^’Calculated from Equation (6) with v given by v from 
Table I times ^Qxpanded Table. 
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Concluaions 


These calculations indicate that the 
difference in power loss between rtodular and 
linear patterned high-voltage arrays is fairly 
small. While the use of modular patterns can 
reduce the effective mean potential by MO 
percent# for the type of configuration consid- 
ered here we also have a 10 percent increase 
in sheath area# leading to only a few percent 
change In total power loss. Correct compar- 
ison of the relative power losses in the 
linear and modular cases requires use of the 
eKpanded resolution code, Comparison based 
on Table I alone incorrectly predicts a re- 
duced power loss for the modular arrangement. 

These preliminary calculations suggest 
that plasma ‘power loss should not bo a primary 
consideration in designing the physical ar- 
rangement of high voltage arrays. 
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PtASMA COLLECTION BY HIGH VOLTAGE SPACECRAFT AT LOW EARTH ORBIT 

1. Katz,* M. J. Mandall,** G. W. Schnuallt,** D. E. Parka, + P. G. Staan^ 

Syitems# Scl«nc« and Software 
La Jollai California 92038 

Abstract plaiinaa Early studies by Kenner ud ^ have 

shown that small pinholes in insulation 


A computer model q£ the three-dimension- 
al sheath formation and plasma current 
collection by high voltage spacecraft has 
been developed. By using new space charge 
density and plasma collection algorithms, 
it is practical to perform calculations 
for large, complex spacecraft. The model 
uses NASCAP compatible objects and geo- 
metries. Results indicate that ion focus- 
ing observed in the laboratory during high 
voltage collection experiments is probably 
due to voltage gradients on the collecting 
surfaces. 

Nomenclature 

0 Renormalized screening length 

e Electron charge 

1 Total collection cur^-ent 

j One-sided electron (ion) thermal 
current 

k Boltzmann constant 

L Object dimension; computational mesh 

constant 

f£^ Variational function for 

m Electron (ion) mass 

n^ Ambient plasma density 

r General space point 

'it Renormalized plasma temperature 

V Particle Velocity 

Particle thermal velocity 

X Distance from sheath edge 

c Functional derivative operator 

Permittivity of free space 

0 Plasma temperature 

Ap Debye (screening) length 

p Space charge density 

Derivative of space charge function 
with respect to potential 
(J) Electrostatic potential 

Solution to Ith linearized Poisson 
equation 

jth iterate of solution for 
I. Introduction 

As satellite electrical power require- 
ments increase, there is constant effort 
to increase operating voltages on solar 
arrays. The prime driver is weight reduc- 
tion in both power cables and elimination 
of DC to DG converters used to supply high 
voltages to high power microwave tubes or 
ion thrusters . One of the primary diffi- 
culties associated with high voltages on 
satellites in low earth orbit is the power 
drain due to leakage currents flowing be- 
tween exposed surfaces through the ambient 

^’Pr^ram Manager , Research Scientist , 
tSenior Research Scientist, ^Junior Re- 
search Scientist* 




can make almost perfectly insulated solar 
arrays couple to an ambient plasma as ef- 
fectively as bare metal plates* Thus just 
covering exposed high voltages may not be 
ah adequate solution. 

In order to predict parasitic losses 
due to surrounding plasma# some simple 
models of plasma collection have been de- 
vised. These include models for current 
collection from other plasmas, such as ion 
engine efflux 2#3 and simple models de- 
signed for LEO. 4,5 Recently large scale 
experiments have been performed in the 
thermal vacuum test chamber at NASA Johnson 
Space Center to simulate collection in LEO 
by a large (1 m x 10 m) solar array. 4 

The purpose of this paper is to intro- 
duce a new three-dimensional computer 
model, NASCAP/LEO, which is designed to 
calculate plasma leakage currents. We 
will present the physical arguments that 
led to the model followed by a brief de- 
scription of the numerical techniques 
employed. We will then show results of 
applying this model to the NAS A/ JSC ex- 
periments of McCoy and Konradi. 4 These 
simulations reproduced the observed ion 
focusing phenomenon. 

II. Theory 

We develop here a simplified nonlinear 
screening model for electric potentials 
and plasma structure surrounding an object 
with exposed high voltage surfaces. We 
define high voltages, pf by requiring 
that the potential energy associated with 
the voltage , ep , be much larger than the 
plasma thermal energy, 0; 

ep >> k0 . 

We also require that the Debye length, X^, 

of the plasma be less than, or at most 
comparable to, typical object dimensions, 

Lj 

Xjj ^ L. 

This clearly places us in the space charge 
limited collection regime. 

In the limit of very small Debye length 
and very high voltages a quasi -one-dimen- 
sional sheath forms in front of the col- 
lecting surface. The current in the 
sheath is the one-sided random thermal 
current: 


jg * n^e(ke/2Trm)^/^ . (1) 



Th« velocity within the sheath is 


( 8 ) 


y mV* ■ j mv| + e(>(r) (2) 

where 

- (k0/2irm)^ '* (3) 

end the char 9 e density is 

P • • ^ • ( 4 ) 

for e^ >> k0 this reduces to 

p ♦ -n_a(k6/4ire(>) (e(> >> k0) (5) 

W • 

When the voltage drops down near the 
plasma potential the charge density is 
determined by the exclusion of oppositly 
charged particles and a Debye screening 
expression is valid: 




4 + p(<^) 


O Ap 


e_k0 

o 


(e<j» < k0) 


( 6 ) 


Laframboise and Parker have shown that 
this expression is valid at very small 
potentials even for collisionless plasmas.^ 
Since the space charge limited potentials 
.4/3 


have a x 


form and the electric field 


has an x^'^^ dependence p then half the 
screening of the field is accomplished in 
the first one-eighth of the sheath. If 
object dimensions are large or comparable 
to sheath dimensions p then the one-dimen- 
iSional charge density formulation can be 
used as an effective nonlinear screening 
to determine sheath potential contours - 

Combining the low and high potential 
expressions we can construct a single 
analytical formula which is a function of 
only the local potential and depends 
parametrically on the neutral plasma tem- 
perature and density: 


PJ±L 

"o 


- + /4r (6(Ji/k0) 


(7) 


This expression has the properties of 
going to the low potential Debye shielding 
limit (equation 6) when e<{> << Jc0 / and to 
equation (5) for e^ >> k9 . The advantages 
of this expression over a particle track- 
ing formulation for obtaining charge den- 
sities are overwhelming. Statistical and 
numerical difficulties preclude the direct 
approach for all but the very simplest 
problems. Using the screening formulation 
presented here reduces the complex Poisson- 
Vlasov sheath problem to a nonlinear equa- 
tion which for small potentials is Helmhols 
in character: 


We then use these potentials surrounding 
an object to define an approximate sheath 
boundary. Typically we define the sheath 
to be the equipotential surface (p » k9/e. 

The sheath current density at the 
sheath boundary is the plasma thermal cur- 
rent. Thus the total current collected 

is 


^ * ^shi^ath boundary ^o 

- j. 


^o " A“*sheath 


( 9 ) 


The distribution of currents on the space- 
craft is found by the forward pushing of 
representative particles from the • sheath 
boundary until they impinge upon the 
satellite. The particles start out with 
thermal velocity normal to the local 
sheath surface and are accelerated by the 
sheath electric fields. By integrating 
the product of the collected current and 
the local surface potential the power 
drain can be computed. 

The nonlinear screening model can be 
made more accurate by modifying the shield- 
ing due to focusing effects. The particle 
tracking algorithms can calculate local 
enhancement of charge density due to the 
convergence of sheath particles. These 
terms can be used to modify the screening 
iteratively to reach a self-consistent 
solution. The final results then would 
be accurate for all ratios of object dimen- 
sions to sheath lengths. 

III. Code Implementation 

The implementation of the above ideas 
has been performed using the basic tech- 
niques developed for the standard NASCAP 
code. Calculations proceed in four 

phases: (A) Object definition; (B) Solu- 

tion of the nonlinear Poisson equation 
(equation 8); (C) Determination of the 
sheath boundary and its area; and (D) 
Particle tracking to determine the cur- 
rent distribution on the test object. 

A. Object Definition 

Objects for NASCAP/LEO are defined 
using the standard NASCAP code. This 
takes advantage of the simple input and 
powerful graphics developed for NASCAP. 
However, the preliminary version of NASCAP/ 
LEO does not handle the full generality of 
NASCAP-defined objects. 

B. Nonlinear Potential Solution 

The nonlinear Poisson equation (equa- 
tion 8) is solved iteratively by succes- 
sive linearization (a multidimensional 
Newton-Raphson method) . The linearized 
form of equation (8) is 
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( 10 «) 

or 

-€o7^i - - 

whoro o' donotos Oo/3^), and donotos 

th« Zth "major** itorato of tha solution. 

Tha linaar aquation (10b) is solvad itara- 
tivaly (laading to a saquanca of "minor" 
itaratas using (as in NASCAP) a finita 

alamant# scalad conjugata gradiant mathod 
with provision for multipla nastad grids 
with compatibla boundary alamants. 

Tha finita alamant formulation for (10b) 
is darivad from tha variational principla 

with 

Two points should ba mads concarning aqua- 
tions (11). First, sinca aquation (7) is 
non-monotonic, o' can ba positiva, laading 
to a non-positiva dafinita sat of aqua- 
tions. Wa hava found it satisfactory to 
circumvant this problam by malting tha ra- 
placemant 

o'(Oj-i) max^O, (11c) 

in (11b) . Altarnativaly , a mathod de- 
signad to handle non-positiva-daf inite 
problans might be adopted. 

Second, if o' is too far negative so 
that the second term of (11b) dominates 
tha first, tha solution of aquations (11) 
is an unphysical oscillatory potential. 

This takes place for < 0.7 L, where L 

is the local mesh spacing. These oscilla- 
tions are eliminated by replacing (7) by 
a formula which has tha same values for 
large but no more screening than tha 
coda can resolve at small 


. . > r, . „/T, 

where 

D ■ max (X^, 0.7 L) 

and 

For the typical case T >> )c0/e, the 


modifiaii screening should ba compensated 
for by defining tha sheath boundary to ba 
at soma potential intermediate between 9 
and T. 

For tha sample calculations discussed 
below, approximately 3-5 minutes of (JNIVAC 
1100/81 time were required for solution to 
aquation (8). 

C. Shaath Boundary Determination 

Tha Shaath boundary is determined by a 
sat of routines which examine each finita 
alamant for tha presence of a user- 
specified Shaath boundary potential con- 
tour. Tha thraa-dimansional sheath con- 
tour is approximated in each selected 
alamant by a collection of contiguous tri- 
angles. Tha entire sat of neighboring 
triangles defines tha sheath boundary to- 
gether with its area, and thus tha net 
currant collected. 

D. Particle Trac)cinq 

Particle trajectories are followed by 
forward trac)(ing from the centroid of each 
sheath boundary triangle. Each particle 
is given an initial currant which is pro- 
portional to tha local ambient flux and 
to the area of the triangle, an initial 
)cinetic energy equal to a user-specified 
value, and an initial direction antiparal- 
lal to the local E-fiald (for electrons). 
The particle trac)cing is performed using 
Boris' second-order leapfrog schema 
with a dynamically adjustable timestap. 
Provisions are included for magnetic 
field affects on particle motions. Tha 
actual particle trac)cing routines of 
NASCAP/LEO are straightforward adaptations 
of the NASCAP emitter particle tracicing 
subroutines . 

IV. Sample Calculations 

To illustrate the capabilities of 
NASCAP/LEO we studied a system similar to 
one used in experiments by .McCoy and 
Konradi. ^ The system consisted of an 
object 8.7 X 1.8 * 0.3 meters with 1.2 < 

1.2 meter patches which can be biased on 
one surface (Fig. 1) . This is intended 



Fig. 1 Test object for sample calcula- 
tions. The labeled portions repre- 
sent a conductive scrip, while the 
border, sides, and back are in- 
sulating plastic. Biasing indi- 
cated is for case (d) . 


(12a) 

(12b) 

(12c) 
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to roprotont « plastic shsst with a con- 
ducting mstal strip. Th« total arsa 

37.6 and ths matal arsa was 10.4 m^. 

It was modslsd with 0.3 m rssolution in a 
volums 9.6 X 9.6 « 19.2 motors, fillod 

with 2.3 oV plasma with donsity 1.3 * 10^ 

cm*^, and thus a Oobyo longth of 1 cm and 
ono-sidod oloctron thormal curront of 

53 mA/m^. 

All simulations woro dono in positivo 
bias (oloctron-collocting) mods. Similar 
rosults would bo obtainod for nogativo 
(ion collocting) objocts, with currants 
roducod by tho squaro root of tho mass 
ratio. Sovon tost casos aro shovm horot 

(a) ontiro objoct at *^100 V. 

(b) ontiro objoct at ^2400 V. 

(c) plastic at ground, motal at >»>2400 V. 
(d* plastic at ground, motal biasod as 

indicatod in Fig. 1. 

(o) samo as (d) , but with potontials 

doublod - 1200 V) . 

max 

(f) samo as (d) , but with potontials 
'quadruplod (V ■ 2400 V) . 

ITlttX 

( 9 ) same as (d) , but with po;.entiais 

OCtupled • 4800 V) . 

max 

Tho shoaths for cnsos (a) and (b) aro 
illustratod in Fig. 2. For tho 100 V case 

tho sheath area was 69 m^ , and to: tho 
2400 V case (b) tho sheath area was 

258 m^, loading to collection currents of 
3.66 and 13.7 amperes respectively. 

In tho actual experiment tho object 
was not at uniform potential, but rather 
tho insulating plastic collected sufficient 
current to remain near plasma ground. This 
was tho motivation for cases (c)-(g), which 
provided more interesting potential con- 
tours and particle trajectories. 

Tho sheath for case (c) is shown in 
Fig. 3. It is apparent that the sheath 
spills over the ends and sides to sub- 
stantially enhance the current collection. 

The kheath area for this case was 92 m^, 
giving a current of 4.9 amperes. Note 
also the electrostatic focusing indicated 
by the x-y cross-section of the particle 
trajectories, tending to concentrate the 
current near the center of the strip. 

The sheaths for cases (d)-(g) are shown 
in Figs. 4-5. The sheath areas are 27, 

37, 61, and 81 m^ respectively. Electro- 
static focusing effects similar to case 
(c) are observed. The current versus 
voltage is shown in Fig. 6 . Except for 
the high-voltage end of each panel, the 
points are not far from a universal curve, 
indicating that a two-dimensional repre- 
sentation would be fairly good in this 
case. Area enhancement factors range 
from \2 at 100 V to 'vlO at 3 kV. The 
greatly enhanced current to t.he end 


T 

I 


100 V 
y -2 view 




2400 V 
y-z view 





Fig. 2 Electron sheath particle trajec- 
tories for cases (a) and (b) , with 
the test object at uniform poten- 
tial . 


section of each panel is a three-dimen- 
sional effect, and makes a substantial 
contribution to the power drain on the 
array. 

A comparison of the results presented 
here with those of McCoy and Konradi must 
be approximate because the plasma in the 
experiment was far from homogeneous. The 
temperature and density we have chosen in 
our calculation correspond roughly to the 
estimated mean experimental values. The 
collection current per unit area from our 
calculation can be converted to Argon ion 
current by dividing by the square root of 
the mass ratio, yielding a one-sided 

plasma current of 0.2 mA/m^. For the 
0-4800 volt bias case, the predicted Argon 
ion current collected per unit area then 

IS 1.6 mA/m^. This is quite close to the 
experimentally observed current of 

2 mA/m^ , far closer than one would expect 
due to the uncertainty in the experi- 
mental plasma parameters. That the cal- 
culation and experiment are in reasonable 
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x-y vimvt 


y-2 view 


x-2 view 


Fig. 3 Electron sheath particle trajec- 
cories for case (c) • (Plastic at 
0 V, metal strip at 2400 V.) 

agreement, however, is certainly encourag- 
ing. 

V. Summary 

We have presented calculations using a 
preliminary version of NASCAP/LEO, a fully 
three-dimensional model for the study of 
current collection by laroe, high-voltage 
spacecraft. These calcula .ions hav^e repro- 
duced both qualitative features such as 
ion focusing and quantitative ion collec- 
tion currents observed in laboratory ex- 
periments. Extensions of NASCAP/LEO to 
include more complex geometries and mater- 
ials as well as finer spatial resolution 
will increase the accuracy and applica- 
bility of the model. NASCAP/LEO should 
be a useful engineering tool for designing 
the spacecraft of the future. 






Fig. 4 y -2 views of electron sheath 

particle trajectories for cases 
(d)-(g). Case (d) is shown twice 
scale. 
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V 

Fig. 6 Current per 1.44 section of 
panel versus potential for the 
four cases (d)-(g). Note sharp 
increase in current to the end 
section, and the relative insensi- 
tivity of section current to other 
sections. The cross indicates the 
mean current to a similar area of 
a panel whose conducting strip is 
uniformly held at 2400 volts 
(case (c) ) . 
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7.3 


FLUX DEFINITION FILE 


There are four basic NASCI&F flux types ~ test tank 
beam^ maxwelXlan> particle pushing; and double maxwellian. 

The first card in the flux file identifies which of these 
types you are using. It must have the word ’T'iPE' followed 
by an integer from 1 to 4. 

This section describes the iijput for each flux type. 
The most complicated case is TYPE X, the test tank electron 
beam. You may skip over that subsection if you don't intend 
to use it. 

Free format input is used throughout NASCAP. The test 
tank beam relative intensity cards for polar coordinates are 
the exception. This one set of cards requires a specific 
format (see below) * All other data can be placed anywhere 
in the first 80 columns of a data card image. 


7.3.1 Test Tank Electron Beam — Type 1 

Type 1 flux is for a test tank with an electron beam 
source. There are actually two input forms, depending on 
whether the user specifies the beam pattern in polar or 
Cartesian coordinates. But these two types of input are 
almost identical; differing in cards number 2, number 10, 
and the pattern specification. 

Card Number Contents 
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1 


TYPE J 


n 


2 

3 

4 


Beam pattern option (see below) 

Position of calibration plane on 
NASCAP grid — constant Z value (1-3 3] 

Beam dimensions (X, Y) in meters 


5 


Beam energy (eV) 


6 


Tracking timestep velocity 
(grid units/timestep) 


7 




Gun location on NASCAP grid 
(X, Y, Z) [outer grid units] 




♦ 



SAMPLB 

LOCATION 


T 


i 


CALIDRATIOH 

plane; 

LOCATION 


Z 


Figure 7.3. Geometry of the test tank configuration. 
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BEAM KIOTH 



Card Number 
8 

9 

10 


-END 


Contents 

Position of sample plane on 
NASCAP grid - constant Z value 
[1-33, inner grid units) 

Total beam current (amperes) 

(Cartesian input only) number 
of points (NX, NY) in calibra- 
tion input grid (2 integers be* 
tween 2 and 33) 

Specify beam pattern: 

Polar — 36 formatted cards. 
Cartesian - any number of 
free format cards NX.NY 


The four things that NASCAP needs to know are (1) Be- 
ginning and end of the beam, (2) Beam intensity and energy, 
(3) Timestep for tracking beam particles, and (4) Beam 
pattern. NASCAP also finds out here if you just want to 
run a beam calibration or want to continue with one or more 
TRILIN iterations. The user does not specify beam direction 
it is assumed to be parallel to the Z axis. 


Beginning and End of Beam 

Cards 7, 3, and 8 are used to specify the beam limits. 
Card 7 gives the electron gun location. The location is given 
in off-center (1-17, 1-33 style) NASCAP grid units, and it 
must be in terms of the outermost grid. 

For the end of the beam you must specify a calibra- 
tion planf^ and a sample plane. The calibration plane is the 
plane at which the beam pattern was measured. The sample 
plane is the point at which particle tracking is discontinued. 
Therefore, the NASCAP object should be located between the 
electron gun and the sample plane. 


^ = ir 


Figure 7.4. 


Specification of electron gun characteristics 
for the polar coordinate test tank case. Elec- 
tron fluxes at the sample plane for an uncharged 
environment are to be specified at the solid 
points. Notice that the Z axis is pointing into 
the figure. We are looking at the calibration 
plane from the point of view of the electron 
gun. 








fc 


I 


I 


For both the calibration plane and sample plane ^ the 
location is specified by a single Z value. This is the plane 
position in inner grid units. 

Beam Intensity and Energy 

Card 9 gives the beam intensity. This is the total 
beam current in amperes. 

Card 5 gives the beam energy in eV. It is a mono- 
energetic beam. 

Timestep for Tracking 

The user supplies the tracking timestep indirectly, by 
specifying how far an electron, having the given beam energy, 
will travel in one timestep. For instance, if the tracking 
velocity is given as .3 grid units/timestep, the timestep will 
be chosen so an electron just leaving the gun will travel .3 
inner grid units before its velocity is recomputed for the 
second timestep. This number goes on card 6. 

Beam Pattern 

The various aspects of beam pattern that the user 
must specify are: choice of polar or Cartesian input, cali- 

bration pattern size, number of calibration points (Cartesian 
input only) , and relative beam intensity matrix (many cards) . 

Card 2 indicates the coordinate system, and tells 
NASCAP whether to stop after calibration or continue with 
the TRILIN cycle. Permitted options are: 

POLAR RUN 

or POLAR CALIBRate 

or CARTES ian RUN 

or CARTES ian CALIBRate. 
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Any first literal other than CARTES will default to POLAR. 

Any second literal other than CALIBR will default to RUN. 

Specification of the CALIBR option will result in 
cleanly exiting NASCAP as soon as test tank particle track- 
ing is completed. 

Card 4 gives the dimensions of the calibration pattern 
in meters in the X and Y directions. If you are using 
Cartesian coordinates, this isi the distance from the first 
to the last grid point in either direction on your rectan- 
gular mesh. If you use polar coordinates, the beam radius 
will be the larger of the two values. 

Card 10 gives the number of calibration points. For 
polar coordinates you must use an 11 x 36 mesh, so this card 
is omitted. For Cartesian coordinates, the lower bound in 
either direction is 2. The upper bound is 33. If your cali- 
bration plane dimensions from card 4 are (X, Y) and your 
number of points is NX by NY, then the distance between 
points in the X direction is X/(NX-1). In the Y direction, 
it is Y/(NY-1). The two numbers on card 10 must be integers — 
no decimal point. 

Cards 10-45 (polar) or 11-END (Cartesian) give the 
relative intensities of the calibration points. These can 
be in arbitrary units, since NASCAP will normalize the 
relative intensities to correspond with the total beam cur- 
rent given on card 9. NASCAP internally stores the beam 
using polar coordinates, so after reading the relative inten- 
sities, either polar or Cartesian, it will print out the 
(R, 9) values of an equivalent polar grid. 

If you are using polar coordinates, you must format 
the relative intensity input (11F4 . 1) . You must supply 
36 cards for 36 angles in 10° increments, and each card 
must have 11 values for the 11 radial positions along a 
single angle.’ The first value on each card corresponds to 
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the exact center of the pattern, so it must be the same on 

th 

every card. A particular (i, j) value will be the i value 
on the j card, and it corresponds to the (R, 8) location 

R » X maximum radius (card 4) 

where the angle 6*0 lies along the positive X axis in a 
right-handed coordinate system. 

If you are using Cartesian ccordinates, the relative 
intensities are free format — like everything else in NASCAP. 
You can have from one to fifteen values on each card. 

The values are read with X as the running index. That 
is, the first NGX values read will correspond to locations 
(1, 1) to (NGX, 1) . You do not have to group the input on 
the cards in any particular way. For clarity, you might 
want to start a new card for each new group of NGX values. 

7.3.2 Maxwellian Probe Approximation — Type 2 

For a maxwellian distribution, the user specifies the 
plasma density for electrons and protons, and the plasma 
temperature for electrons and protons. 

Density is in units of particles per cubic centimeter 
or per cubic meter, depending on whether CGS or MKS units 
are specified. Temperature can be given in KEV, EV, JOULES, 
or KELVIN. 

The first density card will be interpreted as electron 
density. The last density card will be proton density. 
Likewise with temperatures. A card reading "END" ends the 
definition. 
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The cards for a maxwellian flux definition are then, 
the 'TYPE 2* card, one or more density cards, one or more tem- 
perature cards, and an 'END' card (optional). 


7.3.3 Reverse Trajectory Particle Tracking — Type 3 

Card 2 contains one of the two literals MAXWEL or 
DEFOR. Following 'M/^WEL' input is as in Type 2, (q.v.), 
except an END card is required. Following 'DEFOR', cards 
are required giving lUNIT, HOUR, IDAY. The data contained 
in the ELT NASCAP. DEFOR must have been copied onto file 
lUNIT. 


Subsequent cards contain 
NSPEC 
NENG 
NTHET 
NPHI 
NSTP 

RMASS(l), CHARGE (1) 


Number of species 

Number of energies 

Number of polar angles 

Number of azimuthal angles 

Maximum number of steps 
per particle 

Mass (kg) , Charge (C) 


RMASS (NSPEC), CHARGE (NSPEC) 

VCODE Initial code velocity 

STV(3) Vector characterizing 

flux anisotropy 

The present code requires NSPEC = 2, with species 1 
being electrons and species 2 protons. The final card, STV, 
may be omitted, in which case the flux is assumed isotropic. 

7.3.4 Double Maxwellian Flux — Type 4 

Flux definition input file consists of the TYPE card 
specifying type 4, followed by up to four cards specifying 
maxwellian components , and (optionally) an END. Each 
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component card consists of five fields, each field containing 
a literal or a floating-point number. 

The fields of the component cards and their permis- 
sible entries are: 


Field 

Contents 

Permitted Entries 

1 

Species 

' PROTONS ' , ' ELECTRONS ' , ' IONS ' , 
'END' 

2 

Density 


3 

Density Unit 

' CGS ' , ' MKS ' 

4 

Temperature 


5 

Temperature 

Unit 

•EV, 'KEV, 'JOULES', 'KELVIN', 
' ERGS ' 


Note that in the current version of NASCAP the terms 'IONS' 
and 'PROTONS' are synonymous. 
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Example 7.1a. (Flux Definition) . Defines a test tank electron 

beam flux using polar coordinates. 
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Example 7,1b. (Flux Definition) . Defines a test tank 

electron beam flux using Cartesian coordinates. 








1 . 

TYPE 

2 

2* 

1*00 

CGS 

3* 

1 *0 

C6S 


. 10*00 

KEV 


10*0 

KEV 


Example 7.2. (Flux Definition). Defines a neutral 10 keV 
plasma for use in the Maxwell probe operating 
mode . 
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Example 7. 3a. (Flux Definition). Defines a neutral 10 keV 

maxwellian plasma for use in the reverse tra- 
jectory particle pushing mode. 
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Example 7.3b. (Flux Definition). Takes flux from DeForest 

venvironmental data for hour 9.998 of day 73. 
The data is read from LUN 9. 
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Example 7.3c. (Flux Definition). Defines a non-neutral, 

non-equilibrium maxwellian plasma for use by 
the reverse trajectory particle pushing mode. 
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Example 7.4. (Flux Definition). Defines a double maxwellian 
flux environment. 


113 






7,3,5 Multiple Reading of Flux Definition; Changing Environ- 

menta 

NASCAP attempts to read a flux definition each time 
TRILIN or DETECT is called. If the ’DEBYE’ option is speci* 
fled, it is also read when CAPACI is called. If the flux 
file contains only one flux definition, NASCAP will rewind 
the file and re-read the same specifications each time. An 
alternate mode of operation is to construct a flux file con- 
taining several successive flux definitions. This is 
particularly useful for simulation of a changing environment/ 
or for a run having several TRXLIN-DETECT cycles (as DETECT 
requires a TYPE 3 flux specification) . For successive fluxes 
to work properly/ each specification must end Cleanly, so 
that NASCAP encounters a ’TYPE’ card whenever it looks for 
a new environmer»t . Thus types 2 and 4 require an ’END’ card 
if another flux is to follow. Type 3 must have an STV card 
(a zero vector may be used to specify isotropic flux) , or 
an ’@EOF' card. Type 1 must have exactly the correct number 
of data cards. If an end-of-file condition is encountered 
when a 'TYPE' card is expected, NASCAP will rewind the file 
and begin the sequence of environments anew. 

An additional consideration for changing environments 
is that if the ’DEBYE’ option is specified a sudden environ- 
ment change Will cause a sudden potential change due to the 
difference in space charge screening. Also, the output of 
the CAPACI module is environment dependent on the 'DEBYE' 
case. NASCAP accounts for ‘these effects in an approximate 
way to avoid the need to call CAPACI for each new environ- 
ment. The accuracy of the approximations used may be gauged 
by calling CAPACI for the extreme environments used, and 
noting the variation of the spacecraft effective radius. 
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APPENDIX B 


DEFAULT MATERIALS AND PARAMETERS 

DEFAULT MATERIALS 
Conductors 

AQUADG Conducting carbon coating (Aquadag) 

ALUMIN Aluminum 

CPAINT Conducting paint 

GOLD Gold 

INDOX Indium Oxide (In20j) 

MAGNES Magnesium 

SCREEN Metallic screen 

Insulators 

KAPTON Kapton 

NPAINT Nonconducting paint 

SI02 Quartz (Si02) 

SOLAR Solar cell cover material (MgF 2 on SiC 2 ) 

TEFLON Teflon 

Note that the default properties are not necessarily 
the same as those used for the SCATHA model. 
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t0 
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1C 

11 
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13 
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PtJOPCRTY 

crsLECTPic constant 

THICKNESS 

CONDUCTIVITY 

ATOMIC NUHBEH 

DELTA .MAX >COEFP 

t-MAX >OEPTHip*-I 

RANGE 

EXPONENT > RANGE 

range > EXPONENT 
EXPONENT 

YIELD FOR IKEV PROTONS 
MAX DE/DX FOR PROTONS 
PHOTOCURRENT 

surface resistivity 

SPACE DISCHARGE POT’L 
INTERNAL DISCHAR5E?0T*L 


1 

-I 

1 

9 

3 

2 


INPUT VALUE 
l.OC+CO {NONE) 
.CC-C3 
.|CtCD 
.3C*0l 
.7C-01 
.2C*CL 
.6C+02 
1 .3C*0G 
2.4CtC2 
l.73«-CC 
2.4“-CL 
2.30+02 
4.0C-QS 
-l.CC+CQ 
1.0C+C4 
2.0C+C3 
1.70*01 
1.8C*DI 
1.90*01 
2.QC+01 


METERS 

MHO/M 

I NONE 5 

(NONE ) 

KEV 

ANG. 

(NONE) 

ANQ> 

(NONE) 

(NONE) 

KEV 

A/H**2 

OHMS 

VOLTS 

VOLTS 


'MATERIAL 


AQUADG 


MATERIAL 


3; CPA IN T 


PROPERTY 

OIELECTPIC CONSTANT 

THICKNESS 

CONDUCTIVITY 

atomic number 

DELTA MAX >C0EFF 
E-MAX >0EPTH**-1 

RANGE 

EXPONENT > RANGE 

RANGE > EXPONENT 
EXPONENT 

YIELD FOR IKEV PROTONS 
MAX OE/OX FOR PPOTONS 
PHOTOCURRENT 
SURFACE RESISTIVITY 
SPACE DISCHARGE POT’L 
internal OISCHARGEPOT'L 


INPUT VALUE 
3.5C*0a (NONE) 

1. CC-03 .METERS 
-l.CCtCO .'"HO/M 

3.CC*0a (NONE) 
2.1C*CQ (NONE) 
i.SC-Ci KEV 
•l.CC*aa ANG. 

2. CO (NONE) 
1.05* CO ANG. 
9.8C*CQ (NONE) 
3.3C-Oi (NONE) 
l.91*Q2 KEV 
2.CC-05 A/H**2 
•l.QC+CC OHMS 

1.2C+G4 VOLTS 
2.uC*03 VOLTS 
1.7C*C1 
l.8C*0l 
1.90*01 
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CODE value 
l.CC*CC (NONE) 
1.00-02 MESH 
-1.00*00 MHO/M 
1.30*01 (NONE) 
7.21 *00 (NONE) 
1.7 3 -02 ANG -Cl 
3.38*02 ANG. 
4.15*02 ANG. 
l.3D*CD (NONE) 
1.73*00 (NONE) 
2.44-Cl (NONE) 


2.30+02 KEV 
4.00-GS A/M**2 


-3.85-13 V-S/C 
1.00*04 VOLTS 
2.00*03 VOLTS 
1.70*01 
1.80*01 


‘:S 8 : 8 i 



property 

INPUT VALUE 

CODE VAl.UE 

1 

DIELECTRIC CONSTANT 

1.00*00 

(NONE ) 

l.CC+CC 

(NONE) 

2 

THICKNESS 

1 .20-0 3 

METERS 

1. 00-02 

MESH 

3 

CONDUCTIVITY 

-i.cc+no 

MHC/M 

-1.00*00 

MHO/M 

4 

ATOMIC NUMBER 

6.uC*00 

(NONE) 

6.00*00 

(NONE) 

5 

OELTA MAX >COEFF 

l.OG+CC 

(NONE ) 

7. 06 +00 

(NONE ) 

6 

£-?mAX >0£PTHm*-1 

3.00-01 

KEV 

2.21-02 

ANG -Cl 

7 

RANGE 

— 1 * C 0*00 

ANG. 

9 .80 + 02 

ANG. 

3 

EXPONENT > RANGE 

O.QC 

(NONE) 

0.00 

ANG. 

9 

RANGE > EXPONENT 

2.20*00 

ANG. 

1.55+OC 

(NONE) 

10 

EXPONENT 

1.20*01 

(NONE) 

l.uC+OC 

(NONE) 

11 

YIELD FOR IKEV PROTONS 

3 . 30-01 

(NONE ) 

3. 3C-cl 

(None ) 

12 

MAX OE/DX FOR PROTONS 

1.91*02 

KEV 

1.91+02 

KEV 

13 

PHOTOCURRENT 

2.10-05 

A/M**2 

2. 10-05 

A/y*»2 

14 

SURFACE RESISTIVITY 

-1.00*00 

OHMS 

-8.55-15 

V-S/C 

15 

SPACE OISCHAHGE POT’L 

i.CC+04 

VOLTS 

1.G0+C4 

VOLTS 

16 

INTERN A L 0 1 S C H A R 3E ? 0 T • L 

2.00*03 

VOLTS 

2*00+03 

VOLTS 

17 


1 .70*01 


1* 70+01 


13 


1.80*01 


1 .80+01 


19 


1.90*01 


1.90 +01 


20 


2.00*01 


2.00+01 



CODE VALUE 
3.50fCC (NONE) 
1.GC-C2 MESH 
•l.QC*aC MHO/M 
5.CC+CC (NONE) 
3.GS’'G1 (NONE) 
3.41-C2 AN 3 -Cl 
l.CS+03 ANG. 


0. GO 
1.51*00 


ANG. 

(NONE) 


l.QC+CC (NONE) 


3.30-01 (NONE) 
1.91+02 KEV 
2.00-05 A/M**2 
■8.55-13 V-S/0 
1.00*04 VOLTS 
2.00*03 VOLTS 
1.70*01 
1.80*01 
1.90*01 
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MATERIAL 


4; GOLD 


1 

2 

T 

m 

4 

5 

6 
7 
fi 
9 

10 

11 

12 

13 

14 
13 


n 

13 

19 

22 


PROPERTY 

DIELECTRIC CONSTANT 

THICKNESS 

CONDUCTIVITY 

ATOMIC NUM3ER 

DELTA MAX >COEFF 

E-MAX >0EPTH**-1 

RANGE 

EXPONENT > range 

RANGE > EXPONENT 
EXPONENT 

YIELD FOR IKEV PROTONS 
MAX OE/OX FOR PROTONS 
PHOTOCURRENT 
SURFACE RESISTIVITY 
SPACE DISCHARGE ?OT*L 
INTERNAL OISCHAPGEPO T* L 


INPUT VALUE 
l.OO+OC (NONE) 
1.0C-C3 meters 
• l.CO+OQ MHO/M 
7.9C+01 (NONE) 
S.8C-Q1 (NONE) 
a.CC-Ql KEV 
3.30+01 ANG. 
1.63+CC (NONE) 
3.46 + 01 A N G . 
7.QC-01 (NONE) 
4.13-01 (NONE) 
1.3S+Q2 KEV 
2.90-05 A/M**2 
•l.CO+OO OHMS 
1. 20+04 VOLTS 
2. 00 + 03 VOLTS 
I .70 + 01 
1.80+01 
1.90+01 
2.C0+01 


MATERIAL 


INOOX 


4 

T 


4 

B 

6 

7 

5 


11 

12 

13 

14 

}| 

17 

18 


PROPERTY 

DIELECTRI C CONSTANT 
THICKNESS 
CONDUCTIVITY 
ATOMIC NUMBER 
delta max >coeff 

E-MAX >0£PTH+«-l 

RANGE 

EXPONENT > RANGE 

RANGE > EXPONENT 
EXPONENT 

YIELD FOR IKEV PROTONS 
MAX OE/DX FOR PROTONS 
PHOTOCUPRENT 
SURFACE fiESISTlVITY 
SPACE DISCHARGE POT+L 

Internal oischargepot’l 


INPUT VA 
1.2C+0Q 

1 .00-03 
-l.CC+OQ 

2 .44 +01 
1.40+00 

8.30- 01 
•l.CC + OQ 
C.CO 

7 * 13 + 00 
5.S5+01 

3.30- 01 
1.91+02 
3.20-05 


-UCC + CO 


1.00+04 
2.00+03 
1.70+01 
1 .80 + 01 


LUE 

(NONE) 

METERS 

MHO/M 

(NONE) 

(NONE ) 

KEV 

ANG. 

(NONE) 

ANG. 

(NONE ) 

( NONE ) 

KEV 

A/M++2 

OHMS 

VOLTS 

VOLTS 


CODE VALUE 
1. 00+00 (NONE) 

1.00- 02 MESH 
-l.OC+CC MHO/M 

7.90+01 (NONE) 
2.52+CC (NONE) 

1. 90- 02 ANG-Ql 
1.3S+C2 ANG. 
2.42+01 ANG. 

1. 63+00 (NONE) 

7.00- 01 (NONE) 
4.13-01 (NONE) 
1.35+02 KEV 

2.90- 05 A/M**2 
-3.85-13 V-S/Q 

l.CO+04 VOLTS 
2^00+03 VOLTS 
1. 70+01 
1.80+01 
1.90+01 
2.0C+01 


CODE VALUE 
1. 00+00 (NONE) 
1. 00-02 MESH 
•l.OC+OC HHO/M 
2.44+01 (NONE) 
3. 02 +00 (NONE) 
1. 49-02 ANG-Cl 
1.57+02 ANG. 
O.CO ANG. 

2. 01 + 00 (NONE) 
l.oC+CC (NONE) 
3. 30-01 (NONE) 
1.91+02 KEV 
3.20-05 A/M»*2 
•S. 35-13 V-S/0 
l.On+04 VOLTS 
2,00+02 VOLTS 
1.70+01 
1.30+01 


19 


1,90+01 


1 .90 + 01 


2C 


2.0C+01 


2.C0+01 



MATERIAL 6: MA6NES 






PROPERTY 

INPUT VALUE 

CODE VALUE 

1 

DIELECTRIC CONSTANT 

l.OC+CQ 

(NONE) 

l.CO+CC 

(NONE ) 

2 

THICKNESS 

1.C0-C3 

METERS 

1,00-02 

MESH 

7 

CONDUCTIVITY 

-l.CC+GC 

MHO/M 

-I.GO+GC 

MHO/M 

4 

ATOMIC NUMBER 

1 .2C+C1 

(NONE ) 

1 .20 + 01 

(NONE) 

5 

DELTA MAX >CO£FF 

9. 20 -Cl 

(NONE ) 

7 ,02+00 

(NONE ) 

6 

e-MAX >OEPTH»)t*-l 

2.50-01 

KEV 

2.79-02 

ANG-Cl 

7 

RANGE 

-l.OC+GQ 

ANG. 

6.96+02 

ANG. 

S 

EXPONENT > RANGE 

0.00 

(NONE) 

0.00 

ANG • 

9 

RANGE > EXPONENT 

1.74+00 

ANG. 

1 .75+00 

(NONE) 

10 

EXPONENT 

2, 43+01 

( NONE ) 

1 . QC + OC 

(NONE ) 

11 

YIELD FOR IKEV PROTONS 

2, 44-01 

( NONE ) 

2. 44-01 

(NONE ) 

12 

MAX OF/DX FOR PROTONS 

2. 30+02 

KEV 

2.30+02 

KEV 

13 

PHOTOCURRENT 

2,00-05 

A/M**2 

2. 00-05 

A/M»#2 

14 

SURFACE RESISTIVITY 

-1 .00+00 

OHMS 

-5,85-13 

V-S/C 

15 

SPACE DISCHARGE POT* L 

i.QC+C4 

VOLTS 

l.CC+04 

VOLTS 

16 

INTERNAL OISCHARGEPOTVL 

2.00+03 

volts 

2. 00+03 

VOLTS 

17 


1.70+Gl 


1.70+01 


16 


1,30+01 


l.SC+Cl 


19 


1.90+01 


1.90+01 
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11 

12 

13 

14 
13 
16 
17 

15 

19 

20 



M4TEi?lAL 7: SCREEN 

PROPERTY 

DIELECTRIC CONSTANT 

THICKNESS 

CONDUCTIVITY 

ATOMIC NUMBER 

delta MAX >COEFF 

E-MAX >OEPTH*»-l 

RANGE 

EXPONENT > RANGE 

RANGE > EXPONENT 

exponent 

YIELD FOR IKEV PROTONS 
MAX OE/DX FOR PROTONS 
PHOTOCURRENT 
SURFACE RESISTIVITY 
SPACE DISCHARGE POT'L 
INTERNAL DI SCHARGEPOT'L 


MATERIAL 8: KAPTON 

PROPERTY 

DIELECTRIC CONSTANT 

thickness 

CONOUCTIV ITY 
ATOMIC NUMBER 
delta max >C0£FF 
E-MAX > DEPTH#*-! 

RANGE 

EXPONENT > RANGE 

RANGE > EXPONENT 
EXPONEN T 

YIELD FOR IKEV PROTONS 
MAX OE/OX FOR PROTONS 
PHOTO CURRENT 
SURFACE RESISTIVITY 
SPACE DISCHARGE POT’L 
INTERNAL DISCHARGEPOT’L 


MATERIAL 9: NPAINT 

PROPERTY 

DIELECTRIC CONSTANT 

THICKNESS 

CONDUCTIVITY 

ATOMIC NUMBER 

DELTA MAX >COEFF 

E-HAX >0E?TH**-1 

RANGE 

EXPONENT > RANGE 

RANGE > EXPONENT 
EXPONENT 

YIELD FOR IKEV PROTONS 
MAX DE/OX FOR PROTONS 
PHOTOCURRENT 
SURFACE RESISTIVITY 
SPACE DISCHARGE POTVL 
INTERNAL DISCHARGEPOT* L 



INPUT^VALUE CODE VALUE 

1 - n n ^ n n # m niu c \ i n r ^ r n ^ m n 


i.oc^aa 

(NONE ) 

l.OC+CG 

(NONE) 

1 .GC-C 3 

METERS 

1 . DC -02 

MESH 

-l.uC*OQ 

MHO/M 

-l.OO+OG 

MHO/M 

1 .CC+QD 

(NONE ) 

1 . 00 +CC 

(NONE) 

C • u Q 

(NONE ) 

C.OC 

(NONE ) 

i.ac^cQ 

KEV 

1 . 33-01 

ANC--G 1 

l.CC^Cl 

ANQ* 

1 . 50+01 

ANG. 

1 • SO^OQ 

(NONE) 

D • C 0 

ANG. 

D • D C 

AN 6 . 

1 .50 + 0 C 

(NONE) 

l.GC^DC 

(NONE) 

1 .GC+CC 

(NONE) 

C.DC 

( NONE ) 

w . u C 

(NONE ) 

1 • wC'^OD 

KEV 

l.CO+CC 

KEV 

C.CQ 

A/M **2 

3 . CQ 

A/M »*2 

— 1 •wC^Gu 

OHMS 

- 8 . 35-13 

V-S/Q 

1 . 30+04 

VOLTS 

1 .C 0 +C 4 

VOLTS 

2 .D 0 + 03 

VOL'^S 

2 .GQ +03 

VOLTS 

1 . 7 C+C 1 


1 . 70+01 


l.SC+Ol 


1 . 80 +Cl 


1 . 90+01 


1 . 90 +Cl 


2 .CC +01 


2 .C 0 +C 1 



INPUT VALUE 

CODE VALUE 

3. 50+00 

(NONE) 

3.50+QC 

(NONE ) 

1.27-04 

METERS 

1.27-03 

MESH 

l.CG-14 

MHO/M 

1.30-14 

MHO/M 

5.00+00 

( NONE ) 

5.C0+0C 

(NONE) 

2.1C+CQ 

•(NONE ) 

3.CS+C1 

(.NOME ) 

1.5C-01 

KEV 

4.62-02 

ANG-Cl 

-l.DC+OD 

ANG. 

7.73+02 

ANG. 

C.OC 

(NONE) 

C.DC 

ANG. 

1.42+Cu 

ANG. 

1 .51+00 

(NONE) 

9.sc+ac 

(NONE ) 

l.QC+CC 

(NONE) 

3.3C-0i 

( NONE ) 

3.30-Cl 

(NONE ) 

1.91 + 02 

KEV 

1.91+G2 

KEV 

2.GC-05 

A/M**2 

2.Ga-05 

A/M**2 

l.CC+16 

OHMS 

3. as +03 

V-S/Q 

1.CD+C4 

VOLTS 

1.CC+G4 

VOLTS 

2.CC+C3 

VOLTS 

2.00+03 

VOLTS 

1.70+01 


1.70+01 


1.80 + 01 


1.80+01 


1 .90+01 


1.90+01 


2.00+01 


2.C0+C1 



INPUT VALUE 

CODE VALUE 

3.5C+QC 

( NONE ) 

3.50+GC 

(NONE ) 

5. .30-0 5 

METERS 

5 . 30-04 

MESH 

5 .90-14 

MHO/M 

5 .90-14 

MHO/M 

5 . 3 u + 0 u 

(NONE ) 

5 . CC+OG 

(NONE) 

2.1C+0C 

(NONE ) 

3. 35+01 

(NONE ) 

1 .5C-C1 

KEV 

3.41-02 

ANG-Cl 

-l.OC+OQ 

ANG. 

1.CS+C3 

ANG. 

O.CC 

(NONE) 

C.CO 

ANG. 

1, .05 + 03 

ANG . 

1.51+CO 

(NONE) 

9.8C+CQ 

(NONE ) 

l.CC+CC 

(NONE ) 

3.30-01 

(NONE ) 

3.3C-G1 

(NONE) 

1.91+02 

KEV 

1.91+02 

KEV 

2;.GG-05 

A/M**2 

2.CC-C5 

A/M**2 

l.DC+13 

OHMS 

3.55+CC 

V-S/O 

1. 00+04 

VOLTS 

1.00 + 04 

VOLTS 

1.00 + 03 

volTs 

1.GC+C3 

volts 

1 .7C+C1 


1.7C+01 


1.80+01 


l.SC+Cl 


1.9G+G1 


1.90+01 



118 




'-ATE^IAL 1C: SI02 


1 

2 
3 

«f 

3 

3 

7 

S 

9 

IQ 

11 

12 

H 

15 

16 
17 
la 
19 
2C 


P^OPERTY 

OIELECTPIC CONSTAmT 

THICKNESS 

CONDUCTIVITY 

ATOMIC NUMBER 

DELTA max >C0EFF 

e-MAX >0£PTH*«-l 

RANGE 

EXPONENT > RANGE 

RANGE > EXPONENT 
EXPONENT 

YIELD FOR IKEV PROTONS 
MAX OE/OX FOR PROTONS 
PHOTOCURRENT 
SURFACE, RESISTIVITY 
SPACE DISCHARGE POT»L 
internal OISCHARGEPOT’L 


INPUT VALUE 
4.CC+OC (NONE) 
1.27-04 meters 
l.CC-14 MHO/M 

1. CC-^ni (NONE) 
2.4C+CD (NONE) 
4.QC-01 KEV 

•I.EC+CQ ANG. 
D.CC (NONE) 
1.02+CQ ANG. 
Z.CC^Ci (NONE) 
3.3C-C1 (NONE) 
1.91+C2 KEV 

2. GO-C5 A/M*4i2 
l.C0tl9 OHMS 
1. 00+04 VOLTS 


2.CC+03 VOLTS 
1.70+01 
1.8C+C1 
1.90+01 
2 . 00+01 


MATERIAL 11: SOLAR 


1 

2 

3 

4 
3 
6 
7 

5 
9 

IQ 

11 

12 

13 

14 

15 

13 

19 

20 


PROPERTY 

DIELECTRIC CONSTANT 

THICKNESS 

CONDUCTIVITY 

ATOMIC NUMBER 

DELTA MAX >CCEFF 

E-MAX >0EPTH**-l 

range 

EXPONENT > RANGE 

RANGE > EXPONENT 
EXPONENT 

yield for IKEV PROTONS 
MAX OE/OX FOR PROTONS 
PHOTOCURRENT 
SURFACE RESISTIVITY 
SPACE DISCHARGE POT*L 

Internal oischapsepot»l 


INPUT value 
4.QC+0C (NONE) 
1.27-04 METERS 

1. CC-14 MHO/M 
1. 00 + 01 (NONE) 
4.10+00 (NONE) 
4.1C-01 KEV 

vl.OC+OO ANG. 
C.CO (NONE) 
2i. 30 + 00 ANG. 

2. CS+C1 (NONE) 
1.36+00 (NONE) 
4.00+01 KEV 
2.C0-05 A/M**2 

•l.CO+OQ OHMS 

1. QC+C4 VOLTS 

2. CC+C3 VOLTS 

1. TC+Cl 
1.80+01 
1.90+01 

2. CC+01 


MATERIAL 12: TEFLON 


CODE VALUE 
4.C0+CC (NONE) 
1.27-C3 MESH 
1. 00-14 MHO/M 
l.CC+Cl (NONE) 
1.17+01 (NONE) 
7. 11-C3 ANG-Cl 
1.18+03 ANG. 
C.CO ANG. 
1.69+CC (NONE) 
l.CD+OQ (NONE) 
3.30-01 (NONE) 
1.91+02 KEV 
2.0Q-0S A/.M>+*2 
ai. 85+06 V-S/C 

1. uO+04 VOLTS 

2. CD +03 VOLTS 
1.70+Cl 
l.SO+01 
1.90+01 
2.aQ+Ql 


COCE VALUE 
4.C0+C0 (NONE) 
1. 27-03 MESH 
1. CO-14 .MHC/M 
l.CO+Cl (NONE) 
1.96+01 (NONE) 
1.48-02 ANG-Cl 
5.46+02 ANG. 

0. 00 ANG. 

1.69+00 (NONE) 

1. CO+CC (NONE) 
1.36 +00 (NONE) 
4.00+01 KEV 
2.00-05 A/M**2 

•3.35-13 V-S/C 
1. 00 + 04 VOLTS 
2.0G+C3 VOLTS 
1.70+01 
1. 80+01 
1.90+01 

2. CC+C1 



PROPERTY 

INPUT VALUE 

COOE VALUE 

1 

DIELECTRIC CONSTANT 

2.00+00 

( NONE) 

2.00+00 

(NONE ) 

2 

THICKNESS 

1.27-04 

METERS 

1.27-03 

MESH 

3 

CONDUCTIVITY 

l.CC-14 

MHO/M 

1. 00-14 

MHO/M 

4 

ATOMIC NUMBER 

7 . u C + C Q 

( NONE ) 

7,00+CQ 

(NONE 1 

5 

DELTA MAX >COEFF 

3.00 +0 0 

(NONE) 

1. 98+01 

(NONE) 

6 

; E-MAX >D£PTH**-1 

3.0C-C1 

KEV 

2 • u6 — 0 ^ 

ANG-Cl 

7 

RANGE 

3.00+02 

ANG. 

4.77+02 

ANG . 

3 

EXPONENT > range 

1.59+00 

(NONE) 

2.Q0+02 

ANG. 

9 

RANGE > EXPONENT 

1.00+02 

ANG. 

1.S9+CQ 

(NONE) 

lu 

EXPONENT 

2.00+00 

(NONE) 

2.00+00 

(NONE ) 

11 

YIELD FOR IKEV PROTONS 

3. 30-Cl 

(NONE) 

3.30-01 

(NONE ) 

12 

MAX OE/DX FOP PROTONS 

1.91+02 

KEV 

1.91+02 

KEV 

13 

PHOTO CURRENT 

2.00-05 

A/M*«2 

2.C0-05 

A/M>k7(t2 

14 

SURFACE RESISTIVITY 

l.OC+16 

OHMS 

s .85+03 

V-S/Q 

IS 

SPACE DISCHARGE POT’L 

1. 00+04 

VOLTS 

1 .00+04 

VOLTS 

16 

INTERNAL DI SCHARSEPO T» L 

2.C0+G3 

VOLTS 

2.00+03 

VOLTS 

17 


I.7C+G1 


1.70+Cl 


13 


I.SC+CI 


1.80+01 


19 


1.90+01 


1.90+01 





119 




... 


REFEBENCES 


KatZ/ I., J. J. Cassidy, M. J. Mandell, D. E. Parks, 

G. W. Schnuelle, P. R. Stannard and P, G, Steen, 
’’Additional Application of the NASGAP Code. Volume II. 
SEPS, Ion Thruster Neutralization and Electrostatic 
Antenna Models," NASA CR-165350, February 1981. 

S tunnard , P. R., I. iCa t z , M . i7 . Mande 11, G . J . Ca ss xdy , 

D. E. Parks, M. Rotenberg and P. G. Steen, "Analysis 
of the Charging of the SCATHA (P78-2) Satellite," 

NASA CR-165348, December 1980. 

Phys. Rev., B19, p. 121, 1979. 

Katz, I., J. J. Cassidy, M. J. Mandell, G. W. Schnuelle, 

P. G. Steen, D. E. Parks, M. Rotenberg and J. H. Alexander 
"Extension, Validation and Application of the NASCAP Code, 
NASA CR-159595, January 1979. 

Cassidy, J. J. , "NASCAP User's Manual - 1978,” NASA 
CR-159417, August 1978. 

McCoy, J. E. and A. Konradi, "Sheath Effects Observed 
on a 10 meter High Voltage Panel in Simulated Low Earth 
Orbit Plasma," Spacecraft Charging Technology-1978, 
p. 315, NASA Conference Publication 2071, AFGL-TR-79- 
0082, 1979. 

Laframboise, J. G. , "Theory of Spherical and Cylindrical 
Langmuir Probes in a Collisionless, Maxwellian Plasma 
at Rest," UTIAS Report No. 100, June 1966. 


DISTRIBUTION LIST 


National Aeronautics and Space Administration 
Washington^ D. C. 20546 
Attn: W. R. Hudson/Code RP 

D. P. Williams, Ill/Code RS-5 

National Aeronautics and Space Administration 
Ames Research Center 
Moffett Field, CA 94035 
Attn: H. Lum, Jr./M.S. 244-7 

National Aeronautics and Space Administration 
Goddard Space Plight Center 
Greenbelt, MD 20771 
Attn: R. 0. Bartlett/Code 408.0 

A. Kampin sky /Code 727.0 

E. G. Stassinopoulos/Code 601.0 
R, S. Bever/Code 405.0 

Jet Propulsion Laboratory 
4800 Oak Grove Drive 
Pasadena, CA 91103 
Attn: Ray Goldstein 

H. Garrett 
E. V. Pawlik 
Paul Robinson 
i 

National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 
Houston, TX 77058 
Attn: J. E. McCoy /Code SN3 

A. Konradi/Code SN3 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665 
Attn: J. D. DiBattista/M.C. 158B 

J. W. Goslee/M.C. 364 

National Aeronautics and Space Administration 
Lewis Research Center 
21000 Brookpark Road 
Cleveland, OH 44135 

Attn: Head, Mechanics, Fuels & Physics Section/ 

M.S. 501-11 

Technology Utilization Office/M. S. 7-3 
Report Control Office/M. S. 5-5 
Office of Reliability and Quality 
Assurance/M. S. 500-211 
APSC Liaison Office/M.S. 501-3 ' 


1 copy 
1 copy 


1 copy 


1 copy 
1 copy 
1 copy 
1 copy 


1 copy 
1 copy 
1 copy 
1 copy 


1 copy 
1 copy 


1 copy 

2 copies 


1 copy 
1 copy 
1 copy 

1 copy 

2 copies 




DISTRIBUTION LIST (Continued) 


Department of Electrical Engineering 

Pennsylvania State University 

121 Electrical Engineering 

East Building 

University Park, PA 16801 

Attn: J. Robinson 

Department of Physics 

University of California at San Diego 

P. O. Box 109 

La Jolla, CA 92037 

Attn : E . Whipple 

Aerojet Electrosystems Company 
1100 West Hollyvale Street 
Azusa, CA 91720 
Attn: C. Pischer/Dept. 6751 

Aerospace Corporation 
P. 0. Box 92957 
Los Angeles, CA 90009 
Attn: J. R. Stevens 

R. Broussard 
J. F. Pennell 

Beers Associates, Inc. 

Post Office Box 2549 
Reston, VA 22090 
Attn: Dr. Brian Beers 

Boeing Aerospace Company 
P. 0. Box 3999 
Seattle, WA 98124 
Attn: H. Liemohn/M.S. 8G-23 

D . Tingey 

Communications Satellite Corporation 
Comsat Laboratories 
Clarksburg, MD 20734 
Attn: A. Meulenberg, Jr. 

European Space Agency 
ESTEC 

Zwartweg, Noordwijk 
Netherlands 

Attn: John Reddy, P.B. TTM, 1719 8 2883 


1 copy 


1 copy 


1 copy 


1 copy 
1 copy 
1 copy 


1 copy 


1 copy 
1 copy 


1 copy 


1 copy 











DISTRIBUTION LIST (Continued) 

Ford Aerospace and Conununlcatlons Corporation 
Western Development Laboratories Division 
3939 Fabian Way 
Palo Alto, CA 94303 
Attn: D. M. Newell/M. S. G-80 

J. P her son 

GeniJral Dynamics Convair 
Kearny Mesa Plant 
P. 0. Box 80847 
San Diego, CA 92138 
Attn: J. I. Valerio 

General Electric Company 
Valley Forge Space Center 
P. O. Box 8555 
Philadelphia, PA 19101 
Attn; V. Belanger/U-2439 
A. Eagles 

Grumman Aerospace 
Bethpage, NY 11714 
Attn: M. Stauber 

Hughes Aircraft Company 
P. 0. Box 92919 
Los Angeles, CA 90009 
Attn: E. Smith/M. S. A620 

Hughes Research Laboratories 
3011 Malibu Canyon Road 
Malibu, CA 90265 
Attn; Dr. Jay Hyman 

IRT Corporation 
P. 0. Box 80817 
San Diego, CA 92138 
Attn; J. Wilkenfeld 

JAYCOR 

P. 0. Box 85154 
San Diego, CA 92138 
Attn; E. P. Wenaas 

Kaman Science 

1500 Garden of the Gods Road 
Colorado Springs, CO 80907 
Attn; F. Rich 


1 copy 
1 copy 


1 copy 


1 copy 
1 copy 


1 copy 


1 copy 


1 copy 


1 copy 


1 copy 


124 



1 copy 




DISTRIBUTION LIST (Continued) 

Lee W. Parker, Inc. 

252 Lexington Road 
Concord, MA 01742 
Attn : L . Parker 

Lockheed Palo Alto Research Laboratory 
3251 Hanover Street 
Palo Alto, CA 94303 

Attn: i7. B. Reagan/Bldg. 205, Dept. 52-12 

D . P . Cauf finan 

Martin Marietta Corporation 
P. 0. Box 179 
Denver, CO 80201 
Attn: D. E. Hobbs 

K. Killian 

Massachusetts Institute of Technology 

Lincoln Laboratory 

P. 0. Box 73 

Lexington, MA 02173 

Attn: F. G. Walther 

McDonnell Douglas Astronautics Company 
5301 Bolsa Avenue 
Huntington Beach, CA 92647 
Attn; W. P. Olson 

Mission Research Corporation 
1150 Silverado Street 
P. 0. Box 1209 
La Jolla, CA 92038 
Attn; V. van Lint 

RCA Astroelectronics Division 

P. 0. Box 800 

Princeton, NJ 08540 

Attn: H. Strickberger/M.S. 91 

Science Applications, Inc. 

101 continental Building 
Suite 310 

El Segundo, CA 90245 
Attn; D. McPherson 

Science Applications, Inc. 

2860 S. .Circle Drive 
Colorado Springs, CO 80906 
Attn; E. E. O'Donnell 


1 copy 


1 copy 
1 copy 


1 copy 
1 copy 


1 copy 


1 copy 


1 copy 


1 copy 


1 copy 


1 copy 


DISTRIBUTION LIST (Continued) 

Simulation Physics » Inc.* 

41 B Street 
Burlington/ MA 01803 
Attn: R. G. Little 

SRI International 
333 Raven swood Avenue 
Menlo Park/ CA 90425 
Attn: J. Nanevicz 

TRW Systems 
one Space Park 
Redondo Beach/ CA 90278 
Attn: A. Rosen 

G. Inouye 


126 




1 copy 


1 copy 


i copy 
1 copy 


